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Abstract 

Dixon's multipoles for a system of N relativistic positive-energy scalar particles are evaluated 
in the rest-frame instant form of dynamics. The Wigner hyper-planes (intrinsic rest frame of the 
isolated system) turn out to be the natural framework for describing multipole kinematics. Classical 
concepts like the barycentric tensor of inertia turn out to be extensible to special relativity only by 
means of the quadrupolc moments of the isolated system. Two new applications of the multipole 
technique are worked out for systems of interacting particles and fields. In the rest-frame of the 
isolated system of either free or interacting positive energy particles it is possible to define a unique 
world-line which embodies the properties of the most relevant centroids introduced in the literature 
as candidates for the collective motion of the system. This is no longer true, however, in the case of 
open subsystems of the isolated system. While effective mass, 3-momentum and angular momentum 
in the rest frame can be calculated from the definition of the subsystem energy-momentum tensor, 
the definitions of effective center of motion and effective intrinsic spin of the subsystem are not 
unique. Actually, each of the previously considered centroids corresponds to a different world- 
line in the case of open systems. The pole-dipole description of open subsystems is compared to 
their description as effective extended objects. Hopefully, the technique developed here could be 
instrumental for the relativistic treatment of binary star systems in metric gravity. 
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I. INTRODUCTION. 



An important area of research is nowadays the construction of templates for gravitational 
waves emitted by binary systems. Analytically, this can be done within the framework of 
PN approximations by means of essentially non-relativistic multipole expansions for compact 
bodies. On the other hand, since the main emission are supposed to take place in a region 
where the PN approximation fails, it would be desirable to have at disposal a relativistic 
treatment of multipolar expansions as a preliminary kinematical tool for dealing with open 
gravitating systems. This paper focuses just on the construction of a suitable relativistic 
kinematical background by an N-body system as a tool. A preliminary extension of the 
results of this paper to relativistic fluids is given in [1]. 

Our starting point will be the result recently obtained [2] concerning a complete treat- 
ment of the kinematics of the relativistic N-body problem in the rest-frame instant form of 
dynamics [3, 4, 5, 6]. This program required the re-formulation of the theory of isolated 
relativistic systems on arbitrary space- like hyper-surfaces [7]. and has been based essen- 
tially upon Dirac's reformulation [8] of classical field theory (suitably extended to particles) 
on arbitrary space- like hyper-surfaces {equal time and Cauchy surfaces). For each isolated 
system (containing any combination of particles, strings and fields) one obtains in this way a 
re-formulation of the standard theory as a parametrized Minkowski theory [3]. This program 
shows the extra bonus of being naturally prepared for the coupling to gravity in its ADM 
formulation. The price to be paid is that the functions z^{t, a) describing the embedding of 
the space-like hyper-surface in Minkowski space-time become configuration variables in the 
action principle. Since the action is invariant under separate r-reparametrizations and space- 
diffeomorphisms, first class constraints emerge ensuring the independence of the description 
from the choice of the 3-1-1 sphtting. The embedding configuration variables 2;^(r, a) are 
thus the gauge variables associated with this particular kind of general covariance. 

We summarize here the main results that are necessary to understand all the subsequent 
technical developments. First of all, recall that since the intersection of a time-hke world-line 
with a space- like hypcr-surfacc corresponding to a value r of the time parameter is identified 
by 3 numbers a = t]{t) instead of four, in parametrized Minkowski theories each particle 
must have a well defined sign of the energy: therefore the two topologically disjoint branches 
of the mass hyperboloid cannot be described simultaneously as in the standard manifestly 
Lorentz-covariant theory, and there are no more mass-shell constraints. Each particle with 
a definite sign of the energy is described by the canonical coordinates fjilr), Ki{T) while the 
4-position of the particles arc given by (r) = z^{r, rfi^T)). The following 4-momenta (t) 
are Kj-dependent solutions of pj — mf = with the chosen sign of the energy. 

In order to exploit the separate spatial and time reparametrization invariances of 
parametrized Minkowski theories, we can first of all restrict the foliation to space-like hyper- 
planes as leaves. For each configuration of the isolated system with time-like 4-momentum, 
we further restrict to the special leaves defined by hyper-planes orthogonal to the conserved 
system 4-momentum [Wigner hyper-planes). This foliation is fully determined by the con- 
figuration of the isolated system. One gets in this way[3] the definition of the Wigner- 
covariant rest-frame instant form of dynamics for any isolated system whose configurations 
have well defined and finite Poincare generators with time-like total 4-momentum (see Ref . [9] 
for the traditional forms of dynamics). Finally, this formulation casts some light on the 
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long standing problem of defining a relativistic center of mass. As well known, no defi- 
nition of this concept can enjoy all the properties of its non-relativistic counterpart. See 
Refs.[10, 11, 13, 14, 15] for a partial bibliography of all the existing attempts. 

As shown in Appendix A of Ref.[2] only four first class constraints survive in the rest- 
frame instant form on Wigner hyper- planes. The original configuration variables z'^{T,a), 
fji{T) and their conjugate momenta Pfj,{T, a), Ki{T) are reduced to: 

i) a decoupled particle a;^(r), (the only remnant of the space- like hyper-surface) with 
a positive mass = -y/pf determined by the first class constraint — Mgys ~ 0, Mgys being 
the invariant mass of the isolated system. As a gauge fixing to the constraint, the rest-frame 
Lorentz scalar time Tg = is put equal to the mathematical time: — r ^ 0. Here, 
x^(r) is a non-covariant canonical variable for the 4-center of mass. After the elimination 
of Ts and by means of such pair of second class constraints, we are left with a decoupled 
free point {point particle clock) of mass Msys and canonical 3-coordinates Zg — es{xs — ^x°), 

— * 

kg — [16]. The non-covariant canonical x^{t) must not be confused with the 4- vector 
x^(r) = z'^{r, a — 0) identifying the origin of the 3-coordinates a inside the Wigner hyper- 
planes. The world-line a;^(r) is arbitrary because it depends on a;^(0) while its 4- velocity 
i;^(r) depends on the Dirac multipliers associated with the remaining 4 (or 3 after having 
imposed Tg — r ^ 0) first class constraints (see Section II). Correspondingly, this world- 
line may be considered as an arbitrary centroid for the isolated system. The unit time-like 
4- vector u'^'{p,^) = p^/^s (= 5;^(t) after having imposed Tg — r ^ 0) is orthogonal to the 
Wigner hyper-planes and describes their orientation in the chosen inertial frame. 

ii) the particle canonical variables f]i{T), Kiir) inside the Wigner hyper-planes. They are 
Wigner spin-1 3-vectors, hke the coordinates a. They are restricted by the three first class 
constraints (the rest- frame conditions) = X^ili ~ 0. Since the role of the relativistic 
decoupled 4-center of mass is taken by x'^{t) (or, after the gauge fixing — r ^ 0, by an 
external 3-center of mass Zg, defined in terms of x^ and p^ [2]), the rest-frame conditions 
imply that the internal canonical 3-center of mass ^+ = ^com is a gauge variable that can 
be eliminated through a gauge fixing [17]. This amounts in turn to a definite choice of the 
world- line x'^{t) of the centroid. 

All this leads to a doubling of viewpoints and concepts: 

1) The externa/ viewpoint, taken by an arbitrary inertial Lorentz observer, who describes 
the Wigner hyper-planes determined by the time-like configurations of the isolated system. 

A change of inertial observer by means of a Lorentz transformation rotates the Wigner 
hyper-planes and induces a Wigner rotation of the 3-vectors inside each Wigner hyperplane. 
Every such hyperplane inherits an induced internal Euclidean structure while an external 
realization of the Poincare group induces an internal Euclidean action. 

Then, three external concepts of 4-center of mass can be defined by using the external 
realization of the Poincare algebra (each one corresponding to a different 3-location inside 
the Wigner hyper-planes): 

a) the ea;ierna/ non-covariant canonical 4-center of mass (also named 4-center of spin) x'^ 
(with 3-location a), 

b) the ea;terna/ non-covariant non-canonical M0ller 4-center of energy R'^ (with 3-location 
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c) the external covariant non-canonical Fokker-Pryce 4- center of inertia (with 3- 
location ay)- 

Only the canonical non-covariant center of mass x^(r) is relevant to the Hamiltonian 
treatment with Dirac constraints, while only the Fokkcr-Pryce Yf is a 4-vector by con- 
struction. See Ref. [2] for the construction of the 4-centers starting from the corresponding 
3-centers (3-center of spin[14], 3-center of energy [12], 3-center of inertia[13, 14]), which are 
group-theoretically defined in terms of generators of the external Poincare' group. 

2) The internal viewpoint, taken by an observer inside the Wigner hyper-planes. This 
viewpoint is associated to a unfaithful internal realization of the Poincare algebra: the 
total internal 3-momentum of the isolated system vanishes due to the rest-frame conditions. 
The internal energy and angular momentum are the invariant mass Msys and the spin (the 
angular momentum with respect to (t)) of the isolated system, respectively. By means 
of the internal realization of the Poincare algebra we can define three internal 3-centers of 
mass: the internal canonical 3-center of mass (or 3-center of spin) q+, the internal M0ller 
3-center of energy i?+ and the internal Fokker-Pryce 3-center of inertia y+. However, due to 
the rest-frame condition k+ ^ 0, they all coincide: g*+ ~ -R+ ~ y+- As a natural gauge fixing 
to the rest-frame conditions, we can add the vanishing of the internal Lorentz boosts K 
(recall that they are equal to —R^/Mgys)- This is equivalent to locate the internal canonical 
3-center of mass in a — 0, i.e. in the origin x^(t) = z^{t,0). Upon such gauge fixings 
and with — r ~ 0, the world-line x^(r) becomes uniquely determined except for the 
arbitrariness in the choice of x^(0) [u^{ps) — Pg/^s] 

<(r)=x^(0)+w>,)T„ (1.1) 

and, modulo a;^(0), it coincides with the external covariant non-canonical Fokker-Pryce 4- 
center of inertia, x^(r) = x'^{Q) + ^,'^[2]. 

The doubling of concepts, the external non-covariant canonical 4-center of mass 5^(t) 
(or the external 3-center of mass Zg when — r ~ 0) and the internal canonical 3-center 
of mass ^+ ~ replaces the separation of the non-relativistic 3-center of mass due to the 
Abelian translation symmetry. Correspondingly, the non-relativistic conserved 3-momentum 
is replaced by the external ps = ^sks, while, as already said, the internal 3-momentum 
vanishes, k+ ~ 0, as a definition of rest frame. 

In the gauge where = Mgys, Tg = r, the canonical basis Zg, kg, ffi, Ki is restricted by 
the three pairs of second class constraints = X]j=i '^i ~ 0, g*+ ~ 0; so that 6N canonical 
variables describe the N particles like in the non-relativistic case. We still need a canonical 
transformation ifi, Ki i— > 0], 0], pa, tTo [a = 1, .., A^ — 1] identifying a set of relative 

— * 

canonical variables. The final 6N-dimensional canonical basis is Zg, kg, Pa, vr^. To get this 
result we need a highly non-linear canonical transformation [2], which can be obtained by 
exploiting the Gartenhaus-Schwartz singular transformation [18]. 

At the end we obtain the Hamiltonian for the relative motions as a sum of N square roots, 
each one containing a squared mass and a quadratic form in the relative momenta, which 
goes into the non-relativistic Hamiltonian for relative motions in the limit c ^ oo. This 
fact has the following implications: 
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a) if one tries to make the inverse Legendre transformation to find the associated La- 
grangian, it turns out that, due to the presence of square roots, the Lagrangian is a hy- 
pereUiptic function of aheady in the free case. A closed form exists only for N=2, 

mi = 1712 = m: L = —m\J 4 — This exceptional case already shows that the existence 
of the limiting velocity c forbids a linear relation between the spin (center-of-mass angular 
momentum) and the angular velocity. 

b) the N quadratic forms in the relative momenta appearing in the relative Hamiltonian 
cannot be simultaneously diagonalized. In any case the Hamiltonian is a sum of square roots, 
so that concepts like reduced masses, Jacobi normal relative coordinates and tensor of inertia 
cannot be extended to special relativity. As a consequence, for example, a relativistic static 
orientation-shape S0(3) principal bundle approach [19, 20] can be implemented only by 
using non- Jacobi relative coordinates. 

c) the best way of studying rotational kinematics, viz the non-Abelian rotational symme- 
try associated with the conserved internal spin, is based on the canonical spin bases with 
the associated concepts of spin frames and dynamical body frames introduced in Ref . [2] . It 
is remarkable that they can be build as in the non-relativistic case [21] starting from the 
canonical basis Pa, tt^. 

Let us clarify this important point. In the non-relativistic N-body problem it is easy 
to make the separation of the absolute translational motion of the center of mass from the 
relative motions, due to the Abelian nature of the translation symmetry group. This implies 
that the associated Noether constants of motion (the conserved total 3-momentum) are in 
involution, so that the center-of-mass degrees of freedom decouple. Moreover, the fact that 
the non-relativistic kinetic energy of the relative motions is a quadratic form in the relative 
velocities allows the introduction of special sets of relative coordinates, the Jacobi normal 
relative coordinates that diagonalize the quadratic form and correspond to different patterns 
of clustering of the centers of mass of the particles. Each set of Jacobi normal relative 
coordinates organizes the N particles into a hierarchy of clusters, in which each cluster of 
two or more particles has a mass given by an eigenvalue {reduced masses) of the quadratic 
form; Jacobi normal coordinates join the centers of mass of pairs of clusters. 

However, the non-Abelian nature of the rotation symmetry group whose associated 
Noether constants of motion (the conserved total angular momentum) are not in involu- 
tion, prevents the possibility of a global separation of absolute rotations from the relative 
motions, so that there is no global definition of absolute vibrations. Consequently, an iso- 
lated deformable body can undergo rotations by changing its own shape (see the examples 
of the falling cat and of the diver) . It was just to deal with these problems that the theory of 
the orientation-shape S0(3) principal bundle approach[19] has been developed. Its essential 
content is that any static (i.e. velocity-independent) definition of body frame ioi a deformable 
body must be interpreted as a gauge fixing in the context of a SO (3) gauge theory. Both 
the laboratory and the body frame angular velocities, as well as the orientational variables 
of the static body frame, become thereby unobservable gauge variables. This approach is 
associated with a set of point canonical transformations, which allow to define the body 
frame components of relative motions in a velocity-independent way. 

Since in many physical applications (e.g. nuclear physics, rotating stars,...) angular 
velocities are viewed as measurable quantities, one would like to have an alternative for- 
mulation complying with this requirement and possibly generalizable to special relativity. 
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This program has been first accomphshed in a previous paper [21] and then relativistically 
extended (Ref.[2]). Let us summarize, therefore, the main points of our formulation. 

First of all, for > 3, we have constructed (see Ref.[21]) a class of non-point canoni- 
cal transformations which allow to build the already quoted canonical spin bases: they are 
connected to the patterns of the possible clusterings of the spins associated with relative 
motions. The definition of these spin bases is independent of Jacobi normal relative coordi- 
nates, just as the patterns of spin clustering are independent of the patterns of center-of-mass 
Jacobi clustering. We have found two basic frames associated to each spin basis: the spin 
frame and the dynamical body frame. Their construction is guaranteed by the fact that in 
the relative phase space, besides the natural existence of a Hamiltonian symmetry left action 
of S0(3) [22, 23], it is possible to define as many Hamiltonian non-symmetry right actions of 
SO (3) [24] as the possible patterns of spin clustering. While for N=3 the unique canonical 
spin basis coincides with a special class of global cross sections of the trivial orientation- 
shape SO (3) principal bundle, for N > 4 the existing spin bases and dynamical body frames 
turn out to be unrelated to the local cross sections of the static non-trivial orientation-shape 
SO (3) principal bundle, and evolve in a dynamical way dictated by the equations of motion. 
In this new formulation both the orientation variables and the angular velocities become, by 
construction, measurable quantities in each canonical spin basis. 

For each N, every allowed spin basis provides a physically well-defined separation between 
rotational and vibrational degrees of freedom. The non-Abclian nature of the rotational sym- 
metry implies that there is no unique separation of absolute rotations and relative motions. 
The unique body frame of rigid bodies is replaced here by a discrete number of evolving 
dynamical body frames and of spin canonical bases, both of which are grounded in patterns 
of spin couplings, direct analog of the coupling of quantum angular momenta. 

As anticipated at the outset, we want to complete our study of relativistic kinematics for 
the N-body system by first evaluating the rest-frame Dixon multipoles [25, 26] and then by 
analyzing the role of Dixon's multipoles for open subsystems. The basic technical tool will be 
the standard definition of the energy momentum tensor of the N positive-energy free particles 
on the Wigncr hypcrplanc. It will be seen, however, that in order to get a sensible extension 
of this definition to open subsystems, a physically significant convention is required. On the 
whole, it turns out that the Wigner hyperplane is the natural framework for reorganizing a 
lot of kinematics connected with multipoles. Only in this way, moreover, a concept like the 
bary centric tensor of inertia can be introduced in special relativity, specifically by means of 
the quadrupole moments. 

A first application of the formalism is done for an isolated system of A^ positive-energy 
particles with mutual action- at- a- distance interactions. Then the formalism is applied to 
the case of an open n < N particle subsystem of an isolated system consisting of A^ charged 
positive-energy particles (with Grassmann- valued electric charges to regularize the Coulomb 
self-energies) plus the electro- magnetic field [5]. In the rest frame of the isolated system a 
suitable definition of the energy-momentum tensor of the open subsystem allows to define 
its effective mass, 3-momentum and angular momentum. 

Then we evaluate the rest-frame Dixon multipoles of the energy-momentum tensor of the 
open subsystem with respect to various centroids describing possible collective centers of 
motion. Unlike the case of isolated systems, each centroid generates a different world-line 
and there are many candidates for the effective center of motion and the effective intrinsic 
spin. Two centroids (namely the center of energy and Tulczyjew centroid) appear to be 
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preferable due to their properties. The case n = 2 is studied exphcitly. It is also shown that 
the pole-dipole description of the 2-particle cluster can be replaced by a description of the 
cluster as an extended system (its effective spin frame can be evaluated). This can be done, 
however, at the price of introducing an explicit dependence on the action of the external 
electro-magnetic field upon the cluster. By comparing the effective parameters of an open 
cluster of ni + n2 particles to the effective parameters of the two clusters with rii and 77-2 
particles, it turns out in particular that only the effective center of energy appears to be a 
viable center of motion for studying the interactions of open subsystems. 

A review of the rest-frame instant form of dynamics for N scalar free positive-energy 
particles and some new original results on the canonical transformation to the internal 
center of mass and relative variables are given in Section II. 

In Section III we evaluate the energy momentum tensor in the Wigner hyper-planes. 

Dixon's multipoles are introduced in Section IV. A special study of monopole, dipole and 
quadrupole moments is given and the multipolar expansion is defined. 

After the extension of the previous results to isolated systems with mutual action-at-a- 
distance interactions, we study in Section V the behavior of open subsystems of isolated 
systems, the ccntroids which arc good candidates for the description the collective center of 
motion, and discuss the determination of the effective parameters (mass, spin, momentum, 
variables relative to the center of motion) for the open subsystem. 

Some comments about standing problems are given in the Conclusions. 

The non-relativistic N-particle multipolar expansion is given in Appendix A, while in 
Appendix B contains is a review of symmetric trace-free (STF) tensors. The Gartenhaus- 
Schwartz transformation is summarized in Appendix C. Finally, other properties of Dixon's 
multipoles are reported in Appendix D. 
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II. REVIEW OF THE REST-FRAME INSTANT FORM. 



We briefly review the treatment of N free scalar positive-energy particles in the framework 
of parametrized Minkowski theory (see Appendices A and B and Section II of Ref.[2]). Each 
particle is described by a configuration 3- vector r/i(r). The particle world-line is x^{r) — 
z^ir^fjiij-)), where z'^{T,a) are the embedding configuration variables describing the space- 
like hyper-surface E^. 

The foliation is defined by an embedding i? x S — >■ M^, (r, a) 2;''(r, ct), with S an 
abstract 3-surface diffeomorphic to R^. S^- is the Cauchy surface of equal time. The metric 
induced on Et- is SfAsl^;] = z'^rj^i^ZQ, a functional of z^, and the embedding coordinates 
z'^{t, a) are considered as independent fields. We use the notation cr^ = (r, a'^) of Refs.[3, 5]. 
The z'^{a) — dz^{a)/da'^ are fiat cotetrad fields on Minkowski space-time with the Zf''s 

tangent to E,-. The dual tetrad fields are z^{a) — z^i'T, ^) ^si'^^ ^) — ^b- While in 

Ref. [2] we used the metric convention 77^,^ = e(-| ) with e = ±, in this paper we shall 

use e = 1 like in Ref. [3] . 

If we put \/ g{T, a) = —det gAsij-, ^) and -^7(t, a) = —det gfs{T, a), we have 

zI^{t,B) = (^y^p + ^,,7-4)(r,a) (7'^"^^, = 51), nT,a) = (^j=e^^^p^ z? zl zfj{T,a) 

(normal to E^; t^ij.a) = 1; I ^{r , a) {r , a) — 0) and d^z — z!^{t, a) dr d^T^f^ — 
dr [z!^{t, a) 1^{t, a)] ^/-t{T, a) d^a = ^/ g{T, a) drd^a. 

The system is described by the action [3, 4, 5] 



S ^ j dTd^aC{T,a) = J drLir), 



C{t, a) = -Y^ S%a - f]i{T))mi^ grr{r, a) + 2grf{T, a)f]l{T) + gfs{r, B)f]l{T)f]f{T), 

i=l 

^ I ~ 

^('^) = - ^^iy grrir, Viir)) + 2grf{T, f)i{T))f)l{T) + gfs{T, fji{T))f]l{T)f]l{T). (2.1) 



The action is invariant under separate r- and ir-reparametrizations. 
The canonical momenta are 
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P/.(t, a) = -^^^ = J2 ^^^^ ~ Vi{r))mi 

Zr,,{,r, a) + Zf^,{T, a)vil{T) 



N 



drilir) 

grfir, f)i{r)) + Qfsir, f)i{r))r)l{r) 



Vfrr(T, f]i(T)) + 2grf(T, f]i{r))ril (t) + gfs{r, ffi(T))ril{T)rif{r) ' 
{ri(T),K,,(r)}^-5,,5l. (2.2) 

The canonical Hamiltonian He is zero, but there are the primary first class constraints 

N 



T^A^^^) = Pn ^) - ^)^^^{^ ^1 - ^) «^ir {r) His (t) - 

i=l 

TV 

(r, a)Y\r, a) Y.5\B - fj.^T))^^, ~ 0, (2.3) 



so that the Dirac Hamiltonian is — J d^aX^{T, a)7in{T, a), where X'^{T,a) are Dirac 
multipliers. 

The conserved Poincare generators are (the suffix "s" denotes the hypersurface E,-) 
p^^ J d'ap^{r,a), 

Jr = J d'a[z^{T, a)p%T, a) - z%t, a)p^{T, a)]. (2.4) 

After the restriction to space-like hyper-planes, z^ij^a) ~ x'^{t) + b^{T)a^, the Dirac 
Hamiltonian is reduced to Ho — Xfj,{T)H^{T) + \^j,v{t)'H'^^ {t) (only ten Dirac multiphers 
survive) while the remaining ten constraints are given by 
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N N 



i=l i=l 

^ I 

= 5r(r) - mr)K - b';{T)l^]Y.4{T)^mi + /5f(r) - 



i=l 

N 



- mm{r) - h^,m{r)\Y.^l{T)4{r) ^ 0. (2.5) 



Here 5'^'^ is the spin part of the Lorentz generators 

5f = 6^(r) y dVaV'^(r,a) - 6^(r) ^ dVaV'^(r, a). (2.6) 

The condition = that the unit normal = e^ai3yb'?{T)b'^{r)F^{r), P = 1, to the 
hyper-planes be r-indepedent, implies three gauge fixing constraints = const, for the 
first class constraints 'H^^{t) ~ 0. These imply in turn that the 6^(t)'s depend only on 
three Euler angles describing a rotating spatial frame. Therefore fohations with parallel 
hyper-planes have only 7 independent first class constraints. 

When p1 > 0, the restriction to Wigner hyper-planes is done by imposing 6^ = {bi^ = 
bf) ^ b'^ = L'^u=a{Ps) =s), where L^^{ps, =s) is the standard Wigner boost for time-like 

Poincare' orbits. The indices f — r are Wigner spin 1 indices and we have — L>^o{Ps, =s) — 

u>'{Ps) = Pt/^/^s^ bi^: = e(f(ii(p,)) = L^(p„ =,). 

On the Wigner hyperplane [28], z'^{T,a) = x'^{t) + e!^{u{ps)) a^, we have the following 
constraints and Dirac Hamiltonian[3, 5] 
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N N 



i=l 1=1 
AT N 



or 



N 



- 0, Msys = X] V + ' 

AT 

Psys — — ^ ^ 'I'i ~ 0, 
i=l 

Hd = A'^(r)?Y^(r) = A(r)[e,-M,,,]-A(r) 

j=i 



AT 



i'sir) = K(r),i/^} = A.(r){<(r),7i'^(r)}R. 

-A^(t) = -A(T)ii'^(p,) + e^(ii(p,))A,(T), 

xI{t) = A2(t) - A2(t) > 0, Xs- u{p,) = -A(t), 



<(r) =<(0)-m'^(P3) / dn\in) + e'^iuips)) / rfriA,(ri). (2.7) 

7o Jo 

While the Dirac multipher A(r) is determined by the gauge fixing Tg — r ~ 0, the 3 Dirac's 
multipliers A(r) describe the classical zitterbewegung of the centroid a;^(r) = 2;^(r, 0) which 
is the origin of the 3-coordinates on the Wigner hyperplane. Each gauge-fixing x(r) ~ 
to the 3 first class constraints /«+ ~ (defining the internal rest-frame 3-center of mass 
^cm) gives a different determination of the multipliers A(r) [29] . Therefore each gauge- fixing 
identifies a different world- line for the covariant non-canonical centroid x^P'^{t). Of course, 
inside the Wigner hyperplane, three degrees of freedom of the isolated system [30] become 
gauge variables. The natural gauge fixing for eliminating the first class constraints /t+ Ri 
is x(t) = (?com = ^+ ~ [vanishing of the location of the internal canonical 3-center of mass, 
see after Eq.(2.16)]. We have that ^+ imphes \f{T) — 0: in this way the internal 3-center 
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of mass is located in a unique centroid x^f'^^'^(r) — z'^(t, a — 0) [xi^'^^^ — x'^ — u'^(ps)]. 

Note that the constant x^(0) [and, therefore, also :r^(0)] is arbitrary, reflecting the arbi- 
trariness in the absolute location of the origin of the internal coordinates on each hyperplane 
in Minkowski space-time. The centroid x'^{r) corresponds to the unique special relativistic 
center-of- mass-type world-line for isolated systems of Refs. [31, 32], which unifies previous 
proposals of Synge, M0ller and Pryce. 

The only remaining canonical variables describing the Wigner hyperplane in the final 
Dirac brackets are a non-covariant canonical coordinate x^(t) [35] and p^. The point with 
coordinates x'^{t) is the decoupled canonical external ^-center of mass of the isolated system, 
which can be interpreted as a decoupled observer with his parametrized clock {point particle 
clock). Its velocity Xg{T) is parallel to p^, so that it has no classical zitterbewcgung. 

The relation between x^(t) and x^(t) {a is its 3-location on the Wigner hyperplane) is 
[2, 3] 



p,.S:'^ + e,{Sr-S, 



louPsuPs 



) , (2.8) 



After the separation of the relativistic canonical non-covariant external 4-center of mass 
x^(t), the N particles are described on the Wigner hyperplane by the 6N Wigner spin-1 
3- vectors ffiir), ^^{t) restricted by the rest-frame condition Kj^ = ^^j^ ~ 0. 

The various spin tensors and vectors are [3] 



Js^ — ^sP's ^'sPs + ^s"^ -~ ^sP's ^'sPs + "^s"^' 



- [«''(p.)6'^(«(Ps))-«'^(p.)e'^(«(p.))]4^'' + e'^(«(Ps))e'^(«(p.))-S: 

N 



e'^{u{ps))u''{Ps) - e''{u{ps))u'^{ps)\ J2 vl\/mjc^ + + 

i=l 

N 

'e^{u{p,))e:{u{p,))-e';{u{psM{u{p,))] J^vK, 
St"" = e':{u{ps))e^{u{ps))S^\ 



+ 



i=l 



N N 



i=l 



1=1 



1 



QfJ,!/ QfJ,!/ I 



Qi] ^js grs goi 



Ps + \/^s' 



^ ^ N N N-1 

S = S = ffi X Kj ^ fji X fjj - fl+ X K+ = Pa X TTg. 

i+l i=l a=l 



(2.9) 
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Note that while L^'^ = x'^Pg — x'^p'^ and S^'^ are not constants of the motion due to the 
classical zitterbewung (when A(t) 0), both L^'' — x^p^ — x^p^ and 5'^'' are conserved. 

The canonical variables x^, for the external 4-center of mass, can be replaced by the 
canonical pairs [36] 



^^ Ps ' Ps ' I / K 

Zs ^ es{xs-^x"^), ks^—, (2.10) 



P"s e. 



s 



which make explicit the interpretation of it as a point particle clock. The inverse transfor- 
mation is 



5: = \/l + kKTs+'^), 



Xg — h [J- s H j"'S) 



p^^e^Jl + kl 



Ps 



This non-point canonical transformation can be summarized as [e^ — M, 



N -> 



sys 



1=1 



0] 







P^ 















Ki 



(2.11) 
0, = 

(2.12) 



After the addition of the gauge- fixing — r a; [37], the invariant mass Mgys of the 
system, which is also the internal energy of the isolated system, replaces the non-relativistic 
Hamiltonian Hrei for the relative degrees of freedom: this reminds of the frozen Hamilton- 
Jacobi theory, in which the time evolution can be re-introduced by using the energy generator 
of the Poincare group as Hamiltonian 



After the gauge fixings Ts — t ^ [implying A(t) = —1], the final Hamiltonian, the 
embedding of the Wigner hyperplane into Minkowski space-time and the Hamilton equations 
become 



14 



Hd = Msys - A(t) • 



with 

i;^(r) = w'^(p.) + et'(u(p,))A,(r), 



^,(t) = 0. (2.13) 
The particles' world-lines in Minkowski space-time and the associated momenta are 

<(r) = z'^(r,r/,(r))=x^(r) + e(fKp,)K(r), 

pf(r) = Jmf + k^{t)u^{Ps) + e^(u(p.))fi;,,(r) ^ = (2-14) 



The external rest-frame instant form realization of the Poincare generators [40] with non- 
fixed invariants = ~ M"^ —pl^s ~ ~^%s^ > obtained from Eq.(2.9): 



N 

1=1 

^s + V^s+Pi ~t 

sir us ^rsuQu , rir„s,rsucu 



1 + kjzl r ' ^ x°pl - xlJMl^ + m ' 12 .15) 
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On the other hand the internal reahzation of the Poincarc algebra is built inside the 
Wigner hyperplane by using the expression of Sf^ given by Eq.(2.9) [41] 



N 



1=1 

N 

K+ = 0), 
i=l 

^ 1 
J =^fjiXKi, r ^S' ^ ruv guv ^ gr ^ 



2 

1=1 

N 



i=l 

n = -kI^ > 0, 



='2 „ ^2 



= -<M%^ - kI)S, ^ -eMl^S,. (2.16) 

The constraints — Mgyg fti 0, ~ have the following meaning: 

i) the constraint — Msys ~ is the bridge connecting the external and internal realiza- 
tions [42]; 

ii) the constraints k+ ~ 0, together with K Q (i.e. i?+ qj^ {)) [43], imply 
an unfaithful internal realization, in which the only non-zero generators are the conserved 
energy and spin of an isolated system. 

The determination of for the N particle system has been carried out by the group 
theoretical methods of Ref.[44] in Section III of Ref.[2]. Given a realization of the ten 
Poincare generators on the phase space, one can build three 3-position variables in terms 
of them only. For N free scalar rclativistic particles on the Wigner hyperplane with 
Psys = ~ within the internal realization (2.16) they are: 



i) a canonical internal 3-center of mass (or 3-center of spin) 

ii) a non-canonical internal M0ller 3-center of energy 

iii) a non-canonical internal Fokker-Pryce 3-center of inertia y+. 

It can be shown[2] that, due to k+ ~ 0, they all coincide: g+ ~ Rj^ ~ y+. 

Therefore the gauge fixings x(r) = g+ ~ i?+ ~ ^+ ~ imply A(t) ~ and force the 
three internal collective variables to coincide with the origin of the coordinates, which now 
becomes 

= 2.^(0) + u^'{p^)Ts. (2.17) 

As shown in Section IV, the addition of the gauge fixings xi.^) = Q+ ^ R+ ~ ^+ ~ 
also imphes that the Dixon center of mass of an extended object [27] and the Pirani[45] and 
Tulczyjew[34, 46, 47] centroids [48] all simultaneously coincide with the centroid xF"^^''(t). 
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The external realization (2.15) allows to build the analogous external 3-variables Qg, Rg, 

— * 

Yg. Eq.(4.4) of Ref.[2] shows the construction of the associated external 4- variables x^, Yf, 

— * 

i?^ and their locations a, ay, (Tr on the Wigner hyperplane. It appears that the external 
Fokker-Pryce non-canonical covariant 4-center of inertia F/* coincides with the centroid 
(2.17). 

In Ref . [2] there is the definition of the following sequence of canonical transformations 



v+ 


Pa 







Q+ 


Pqa 


K+ 


T^qa 



1,..,7V-1, (2.18) 



leading to the canonical separation of the internal 3-center of mass (g+, «:+) from the internal 
relative variables pqa, T^qa- Since the rest-frame condition k+ pa implies [2] pqa ~ Pa, 
TTqa ~ in the gauge ~ and in terms of the associated Dirac brackets we get an 
internal reduced phase space whose canonical basis is Pqa = Pa, T^qa = T^a, — ^, ^ — 1. 

The intermediate linear point canonical transformation in (2.18) is [actually this is a 
family of canonical transformations, since the 7oj's are any set of numbers satisfying — 

0) So T'"' ^'^j ~ ~ IV' Taj 7&i ~ ^ab] 



ViV — 



V+ = 



i i 



(2.19) 



The second canonical transformation has been defined in Section V of Ref. [2] by using 
a singular Gartenhaus-Schwartz transformation (see Appendix C), but it was not written 
explicitly for k+ 7^ 0. By using Eqs. (C14) and (C15) we get the following results: 

1) For N = 2 (711 = —712 = we have 
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2' ' 2 

77+ = ^ (771 + 772), K+^Ki + K2, 

P= 771-772, 7f=^(Ki-K2), 

J = r/i X + 7/2 X = 17+ X + 5 = g+ X + S = pX TT, Sq = Pq X TTq, 

^Jm\ + Rl + ^ml + Rl 2 ^^77^2 + ^2 _^ _^ -2 

g+ = i?+ + 

+ ^ ^ ; 

(2.20) 

Then after some straightforward algebra we get (note that in Ref . [2] we used the notation 

M^free) = M^ys and A^(/^ee) = ^) 
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A^(/ree)' 



,2 _L _ . _L r?2 _ ^ ^9 ' ^+ _L ^1 "^2 /^2 



+ — -y / m| + 



+ 



(free) ^ 



+ (mi - ma) W>1 



(free) 



= TTq + F K+, 



P^Pq 



K+ ■ Pg TTg 



B 



Mifree) \l M^^^^^ + K', 



■"^^ Pq + Clfq, 



M(free) ■ 



M 



{free) 



TTg = TT 



K+ 



■-=2— W'-ee) 

A, I 



\def^ 7-1 ^ 

\ — t: — D 



Pq^P + 



AR^ ■ p 



Ml 



(free) 



TT, 



Sq^ S - D pX K+, 



q+ ^ R+ + G Sq X 



Mifree) (Mifree) + Mifree) " 4) ^-^(/.ee) + 4 ^Ifree) + 4 + A^(/.ee)) ' 

Pq + CTTq - G Sq X K+, 



E 



5 ((-)■" - 



A1 



(free) 



+ Kt 



.) (p. + Ctt,) -G-S^, x«+, 



AC,; 



(2.21) 
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2) For N > 2 the results concerning the coordinates (and also g+) are much more involved 
due to the complexity of Eqs. (C15) so that we give only the following results for the 
momenta 



M, 



(free) = \/ + = ^•^(/r-ee) + ^ ~ M(free) = XI J"^? + ^ XI 7ai ^bi T^qa ' TT^b, 

ah 



M{}ree) ^ ^ 



+ 



jmj + iV X Ibi ^qa " ^qh ^A<(/.ee) + 4) > 

y a& 



— * 

«i = + ViV 2J 7ai [TTqa + 



— 1 ) =175 + 



-»2 



M 



{free) 



Y lai Jmf + N "Y ^ai lU T^qa ' T^qb 



ab 



(2.22) 



20 



III. THE ENERGY-MOMENTUM TENSOR ON THE WIGNER HYPERPLANE. 



A. The Euler-Lagrange Equations and the Energy-Momentum Tensor of 
Parametrized Minkowski Theories. 

The Eulcr- Lagrange equations associated with the Lagrangian (2.1) are (the symbol 
means evaluated on the solutions of the equations of motion) 



/ 8r 8r \ 



dL _ dL _ 

where we have introduced the energy-momentum tensor [here ?)j^(t) = (l;f/j(r))] 

T^'M = = -E^'(--^^(-))^====^(-'-)- (3-2) 

Because of the delta functions, the Euler-Lagrange equations for the fields z^{T,a) are 
trivial (0 = 0) everywhere except at the positions of the particles. They may be rewritten in 
a form valid for every isolated system 

a^r^^4 = - ^a^[v^4]r^^. (3.3) 

When ^A[^/gz'^] = 0, as it happens on the Wigner hyper-planes in the gauge ~ and 
Ts — T f=:i 0, we get the conservation of the energy-momentum tensor T"^^, i.e. 9^7"^^ = 0. 
Otherwise, there is a compensation coming from the dynamics of the surface. 

On the Wigner hyperplane, where we have 



xf(r) = z^{T,UT))=x>:{T)+e^,{u{p,M{T), 

xl{r) = grrir, r/i(r)) + 2grr{r, ffi{T))fjl {t) + ^„(r, m{T))f]l{T)f}l{T) = 
xI{t) + 2xs,{r)e'^{u{psM{r) - ff.ir), 

Piir) = ^m2-7-(r,r7,(r))/^i,(T)«:i,(T)P(r,r7i(T))-«:,,(r)7'-^(r,77,(T))^,^(T,77,(T)) = 
= ymH^fW^*''(pJ + er(«(pJX(^) ^P^^rnl 

r ^ 

= d?ap,{T,a)^Y.J^,{T), (3.4) 
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the energy-momentum tensor T^^{t, a) takes the form 



N 



1=1 



i=l 

N 



rrii 






mirjlir) 




^ xI{t) + 2xs^{T)e^r{u{ps)) 








J x^si^) + 2x,i,{T)ef^{u{ps)) 





B. The Energy-Momentum Tensor in the Standard Lorentz-Covariant Theory. 

With the position x^ — z^{T,a), the same form is obtained from the energy momen- 
tum tensor of the standard manifestly Lorentz covariant theory with Lagrangian Ss — 
J drLsir) — — '^f^i'iTii J dr-^/^irj, restricted to positive energies[51] 

-•—1 V "^i \'^^) 



i=l 

-•AB/ 



= e'X{uiPs))eUuiPs))T^''iT,a). (3.6) 
1) On arbitrary space-like hyper-surfaces we have 
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AT 



E 



[^^^(ri,r/i(ri)) + 2;,^(ri, rf,(ri))77[(ri)] 
[<(Ti,r7;(ri)) + <(Ti,77i(ri))r)[(ri)] = 

i=l I'^'lJ 

[z,^(ri,rf,(rO)<(ri,r/,(ri)) + 

+ <(Tl,?7i(Tl)X(Ti,?7i(Ti))?7[(Ti)77,'(Ti) = 

[;^,^(T,r;,(Ti))<(r,r;,(r)) + (<(r, r;,(r))<(r, //.(r)) + 



E 



1 \/a;f(r)v/^ 



r. (T 



[;2,^(r, r;,(ri)) <(r, r;,(r)) + (;.,^(r, r;,(r))<(r, r;,(r)) + 
+ <(r, r/,(r))<(r, 77,(r)))7)[(r) + <(r, 77,(r))<(r, mmrmr)] 

since ciei |2r^| = = ^/ liQrr - Y^grrgrs), 7 = Metgr^sl. 
2) On arbitrary space-like hyper-planes, where it holds[3] 
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9rs — ^rsj T 
= «/rr + X] 9rr 



Qtt 



7 = 1, 



ds^ = xl{T)dT^ + 2i;,(T) • er{u{ps))dTda'' - da^, (3.10) 



we have 




mj(53(a - f)i{T) 



[{x>:{t) + b'f{TX{T)){x:{T) + K{T)vt{r)) + 

+ Kmmmir)]. (3.9) 



3) Finally, on Wigner hyper-planes, where it holds[3] 



we have 
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AT 



rrii 



N 

E 



is(r) . u{ps) ^i2^r) + 2xsp{T)e?{u{psM{ 

+ (i^(r)6,^(«(p.)) +i:r(r)e^f(«(p,)))53(a-r/,(r))r)[(r) + 
+ 6^f(«(p,))6:(«(p,))53fa-rf,(r))7)[(r)r),^(r)l. (3.11) 



Since the volume element on the Wigner hyperplane is u^{ps)d^u^ wc obtain the following 
total 4-momentum and total mass of the N free particle system (Eqs.(2.7) are used) 



P!^ = j dVy-lxf (r) + («(p,))<T»]«„(p.) 



= E 



N 



[ - XiT)u^iPs) + mr) + yir)]e>;^iu(p,))\ 

N N 

\xir)=~i E [V^fc^ + i^Kr)u^iPs) + <(r)6(f(x.(p,))] = J^P'^^^^P^^ 

i=l i=l 

N 



Msys = P^U^iPs) = -A(t) E 



2 



Wr)=-i f;VmH^(r), (3.12) 



i=l 



which turn out to be in the correct form only if A(t) = — 1. This shows that the agreement 
with parametrized Minkowski theories on arbitrary space-like hyper- surfaces is obtained only 
on Wigner hyper-planes in the gauge — r ~ 0, which indeed implies A(r) = —1. 
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C. The Phase-Space Version of the Standard Energy-Momentum Tensor. 

The same result may be obtained by first rewriting the standard energy-momentum 
tensor in phase space and then imposing the restriction pf (r) = 

y /mf - Y'iT, 'niiT))KiriT)Kis{T)l^'{T, f)i{T)) - Kir{T)Y' {t, f)i{T))z^{T, f)i{T)) ^ 

•\/m? + Kf{T)u^^{ps) + Kle'j^{u(ps)) [52] [the last equality refers to the Wigner hyper- 
plane, see Eq.(2.14)]. We have: 



^ 1 r I 

T^^{z{t,B)) = E - / ^riJi:f(rOpr(riK(Ti)<^'(^i(ri)-^(r,a)) = 
on Wigner' s hyper — planes 



E'^'(^-^"'^(^)) 



N 



+ 2xspei{u{ps))r]f - 



r) 



rriiyjxs- u{ps) 



+ 



AT 



E^'(^-^"'^(^)) 



i=l 



vrii L V 

+ A;[(T)(K'^(p,)e,^(K(p,))+i.^(p,)e(f(i.(p,))) + 

<(r)/v?(r) 



+ 



^Jm\ + ref (r) 



:e^f(«(p.))e:(«(p,)) 



(3.13) 



The total 4-momentum and total mass are then 
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1=1 



Xg ■ u{ps) rrii 



Msys 



/\\ I j lib" 



Since in this case we have m^/ ^1 — [f/j(r) + A(r)]^ = Y^rr7|c^~+~^|(r) , the equations above 
show that the total 4-momcntum evaluated from the energy-momentum tensor of the stan- 
dard theory restricted to positive energy particles is consistent [53] with the description on 
the Wigner hyper-planes with its gauge freedom A(r), A(r), provided one works with Dirac 
brackets of the gauge Tg = r, where one has A(t) = —1 and 



x^{t) 



x':{0) + ru^{Ps) + e';:{u{ps)) f dnXr{ri), (3.15) 

Jo 

Therefore, for every A(r), we have 
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N 



+ 



1=1 



<(t;)/v;(7; 



+ 



e^,{u{p,M{u{Ps)) 



T^^T,,a) = Y.5\a-UTs))^rni + Kl{T,)l 

i=l 

N 



N 



i=l 



<(T,)<(r,) 



TV 



•4 = 1 



At/ 



i=l 



(3.16) 
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IV. DIXON'S MULTIPOLES FOR FREE PARTICLES ON THE WIGNER HY- 
PERPLANE. 

In this Section we shall define the special relativistic Dixon multipoles on the Wigner 
hyperplane with — r = for the N-body problem [54] [see Eqs.(2.13) with x^(t) = 

_|_ u'^{j)^) Tg + e'^{u{ps)) jl dTi Ar(ri) = x^s^^^{t) + JJ" dri \{ti)\. By comparison, a list of 
the non-relativistic multipoles for N free particles is given in Appendix A . 

Consider an arbitrary time-like world-line w^{r) — z^{t, ff{T)) — x'^{T)+e!^{u{ps)) rfir) — 
x^s*^^{t) + e'^{u{ps)) ffir) [?7^(t) = ?7^(t) -|- JJ" dri Ar(Ti)] and evaluate the Dixon multipoles 
[25] [61] on the Wigner hyper-planes in the natural gauge with respect to the given world-line. 
A generic point will be parametrized by 



(4.1) 



'-Jo 

so that SZf^{T, a)u^{ps) — 0. 

While for 77(r) = [r/(r) = dTiXr{Ti)] we get the multipoles relative to the centroid 
x^(r), for ff{T) = we get those relative to the centroid x^f^^^{T). In the gauge Rj^. ^ q+ ^ 
y+ ~ 0, where A(r) = 0, it follows that ?7(r) = fj{T) — identifies the barycentric multipoles 
with respect to the centroid xF"^^^(t), which now carries the internal 3-center of mass. 



A. Dixon's Multipoles. 

Lorentz covariant Dixon's multipoles and their Wigner covariant counterparts on the 
Wigner hyper-planes are then defined as 
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{lJ,l...lJ,n){flv) 



{Ts,ff) 



N 



i=l 



+ e(f(i.(p,))e:(i.(p,)) 



N 



+ K(p.)er(«(Ps))+«^(Ps)e^(«(Ps))] 

i=l 



^ / 



i=l 



^ ^/mf + Ki{Ts) 



1=1 



N 



i=l 



= e^^l{u{ps))...e>:::Hps)) 



2 



i=l 



= e^f;(«(p,))...e^:(«(p,)) 



Related multipoles are 



= e>^l(u(ps))....e'^:(u{ps)) 

N 
i=l 

uMpf^^-''"''{Ts,ff) = 0, 

pti}-'^"''{Ts,rf)uM=it^'-'"{Ts,ff), 

n^O ^ p^(r„ ff) = 6^(«(p,))g^-(T,) = P^^ p^. (4.3) 
The inverse formulas, giving the multipolar expansion, are 



n=0 



ni 

n=0 



an 

S^a-rfiT,)). (4.4) 



Note however that, as pointed out by Dixon [25], the distributional equation (4.4) is valid 
only if analytic test functions are used, defined on the support of the energy-momentum 
tensor. 

The quantities q^'-''-^^{Ts,f]), q'^'-''-''^ {T„ff) = g^.^ ' (T„ r/), g;i-™(T„r/) are the 
mass density, momentum density and stress tensor multipoles with respect to the world- 
line w^{Ts) (barycentric for 77 = 77 = 0). 
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B. Monopoles. 



The monopoles correspond to n = [64] and have the following expression [65] (see 
Appendix C for the definition of -^a-^oo) 



Qta 



N 

f^U V 



N 



i=l y de 



^de ^dilei^qd ' '^qe 



N , 1..N-1 N 



(1t{Ts, rf) ^c-»oo XI "^i^^ + ^ XI ^^'"^'" T^ga " T^qb + 0{l/c) = 

«=1 ao 1=1 

AT 

= Y^T^if? + Hrel,nr + 0{l/c), 
i=l 

qJ^{Tsi if) — ~ 0, rest — frame condition {also at the non — relativistic level), 

1..N-1 N 
ah 1=1 

l..Ar-l I..N-1 

= E K^<aK^ + o{ilc) = E + ^(Vc), 



ab ab 



(T3,r;) = t^^(T„r;) = E 



1 s/ml -\- K, 



•2 



N 9 



Y 

1^1 Ede IdileiT^qd ' T^qe 

N 1..N-1 N 

2_^miC -- ^ 2^ 7rg„-7rgb + C'(l/c) 



2 ^ ^ mi 

1=1 ao 1=1 

AT 

= E^^c'-^--'."- + ^(Vc)- (4.5) 

i=l 

where we have exploited Eqs. (5.10), (5.11) of Ref.[2] to obtain the expression in terms of 
the internal relative variables. 

Therefore, independently of the choice of the world-line w^{t), in the rest-frame instant 
form the mass monopole is the invariant mass M = J2iLi V^i ~^ '^h while the momen- 
tum monopole q^ vanishes and q^^ is the stress tensor monopole. 
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C. Dipoles. 



The mass, momentum and stress tensor dipoles correspond to n = 1 [66] 



N 



N 



+ 



The vanishing of the mass dipole q^'^ imphes f]{T) — 77 (r) — JJ^o?TiA(ti) = and 
identifies the world-hne 'w^{t) — x'f'^^'^{T) + e^{u{ps)) -R+ + /J^ dri Ar.(Ti)j . In the gauge 

R+ ^ q+ ^ y+ ^ 0, where A(r) = 0, this is the world-hne w^(r) — xi^'^^'^(T) of the centroid 
associated with the internal M0ller 3-center of energy and, as a consequence of the rest 
frame condition, also with the rest-frame internal 3-center of mass q+. Therefore we have 
the implications following from the vanishing of the barycentric (i.e. A(r) = 0) mass dipole 



g5r^(7;,77) = el\{u{p,))iii,'{T,,fi} = M [rI{T,) - 77^(7;)] = 0, and A(t) = 0, 

^ 77(7;) = ^(r,) ^R+^ q+ ^ y+. (4.7) 



In the gauge R+ ~ ~ y+ ~ 0, Eq.(4.7) with ff = fj = implies the vanishing of the 
time derivative of the barycentric mass dipole: this identifies the center- of-mass momentum- 
velocity relation (or constitutive equation) for the system 

dq^^^{T„rf) ^ _ • ^ ^ 
dTs 

The expression of the barycentric dipoles in terms of the internal relative variables, when 
fj — 7] = Rj^ a; and ^ 0, is obtained by using the results of Appendix C. 
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0, 



qr{Ts,R+) 



N 



N-1 



a=l 



i=l 

Af-1 

\^ r « 

a=l 
N~l 



l..A?"-l 



a=l 

N 



ab 



N 



i=l 



N N-1 



i=l 



N 



1=1 a=l 1=1 



AT-l 



l..iV 



El..iV-l u V 



Efc=l V^^fc + ^ Ede IdklekT^qd 



/fc= 

L.iV iV-1 

EE 

I. .N-1 



^qd ■ TTqe^ 



qa 



'iai iaj r^^i^ 



1..7V-1 



he 



E [^E 

o6c 1=1 



TV v^TV 

laCibilci Z^j=l"^i7ai 



AT 



8=1 

2 



mi 



P>Xc + 0(l/c), 



i=l 
N-1 



q^/A{Ts,R+) = J2(^l-Rl)^ 



I..N 



X 



L.ATTV-l 



Efe=l + ^ Ede Idklek^qd ■ 7?, 



^^y]V(7„i -7„^)p^^^ lbilciT^qh-f^qc + 0{l/c) 

ij a=l 



qe 

1..N-1 



L.AT-l Af 

= ^ E [^E 

abc i=l 



be 
N 



lailbilci Ej=l f^jlaj ] 



m 



P^vr,,-7f,e + 0(l/cX.4.9) 



The antisymmetric part of the related dipole x^^iTgi ff) identifies the spin tensor. Indeed, 
the spin dipole [67] is 
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= M[R:,{Ts) - v'{Ts)] [e>^{u{ps)X{ps) - e';{u{ps))u'^{ps) 

N 



+ 



i=l 



U„{Ps)ST{Ts)[f)\ = 0, 



?7 = R+, 



JJ. bary centric spin for f) — f] — 0, seeEq{2.9), 



1 \/l-fk{Ts 



+ 



+ E 



1^, 



-1 \^l-fk{Ts) 



i 

+ e^^^S:e>^(u(p,))e':(u(p,)). (4.10) 



This explains why "^■u(p^)(7s, ^) is also called the mass dipole moment. 

We find, therefore, that in the gauge /?+ ^ q+ ^ y+ ^ with = Mu^[ps) = 
M x^f'^^^{Ts) the M0ller and barycentric centroid x^f^'^^{Ts) is simultaneously the Tulczyjew 
centroid[34:, 46, 47, 71] (defined by S'^^ Pi, = 0) and also the Pirani centroid[4:5, 72] (defined 
by S^^ xfv^'' = 0). In general, lacking a relation between 4- momentum and 4- velocity, they 
are different centroids [73] . 

Note that non-covariant centroids could also be connected with the non-covariant external 
center of mass and the non-covariant external M0ller center of energy. 



D. Quadrupoles and the Barycentric Tensor of Inertia. 

The quadrupoles correspond to n = 2 [75] 
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1=1 

iV 



N 



and when the mass dipole vanishes, i.e. 77 = Rj^ = fji \/rri^^^r^ / M , we get 



i=l 

AT 



i=l 
AT 



^ l..Ari..7V-l 

= E (7a. -7a,). (4.11) 

ijk ab 

Following the non-relativistic pattern, Dixon starts from the mass quadrupole 

^ I 

?---(T„ R^) = Ek'^lV^? + N'](r.) - mr:^ (4.12) 

i=l 

and defines the following barycentric tensor of inertia 
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u 

^ I 

i=l 

I..N-1 ^ I..N 



ah ijk 
(EiLl + ^ Ede IdhlehT^qd ' TT^e)^ 

de y de 



l..iV-l l..iV 
ab ijk 

^ N ^ 1 / ^Lcd 7ci7di7^qc ■ T^qd ^ ^ Led 7cj7rfi7rgc " TT^^d ^ 



^ Ecd IckldkT^qc ■ T^qd _ ^ l^cd IchldhT^qc " T^qd \ ^ /^2n1 



^2 

/i=i 

L.AT-l 



= n''[qnr]+0{l/c). (4. 

Note that in the non-relativistic hmit we recover the tensor of inertia of Eqs.(All). 
On the other hand, Thome's definition of barycentric tensor of inertia[76] is 
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u 

N 



1..N-1 1..N 



ab ijk 

[Pqa ■ PqbS^^'^' " P^aPll] 
I..N-1 I..N 

Yl Yl ^ (^«^ - laj){lbi - lbk)Pqa ' Pqb^'''''^ " P^^Pg^] X 

ab ijk 

^ l^cd IckldkT^qc ■ T^qd _ j_ ^ Z^cd 1 chl dhT^ qc ' qd \ ^ /(,2n1 



2mf 2mj 



^2 

h=l 

I..N-1 



= E ^a^P^a • P,*>5^^" - P;>;^.] + 0(1/C) = 

ab 

= /'•^'•n?n.]+0(l/c). (4. 

In this case too we recover the tensor of inertia of Eq.(All). 

Note that the Dixon and Thorne barycentric tensors of inertia differ at the post- 
Newtonian level 



^ l..Af-l 1..N 

1 \ ^ \ ^ rrijmk 



thorne 



Y Y ~ laJ){lb^ - Ibk) 



c — ' — ' Nm^ 

ab ijk 



\Pqa ■ PqbO - PqaPqb 



Z^cd laldiT^qc ■ T^qd 



rrii 



+ 0{l/c'). 



E. The Multipolar Expansion. 

By further using the types of Dixon's multipoles analyzed in Appendix D as well as the 
consequences of Hamilton equations for an isolated system (equivalent to d^T^^" = 0) , it turns 
out that the multipolar expansion (4.4) can be rearranged [see Eqs.(Dll)] in the following 
form 
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+ E ^/^^-^"^-(r., rf)6;^^ (^P.))..6;: {u{ps)) Q^,^ Q^,^ S^a - TO), (4.15) 



n=2 



where forn > 2 and ^/ = I^'-''"'"' (T,) = MIl_llJ(i^^■■^^«-^M^^">(^Ts), with J^^-^"^"^"(T,) 
being the Dixon 2^+"-pci/e inertial moment tensors given in Eqs.(DlO). With this form of 
the multipolar expansion, the quadrupole term (n = 2) has the form [see Eq.(Dll)] 



+ <M(Ps))CMPs)) [5?™(T.,r/) - \ (gf^^^"^)"^(T„r;) + 4™)-(T„r/)) + 
+ gf™-)(r„r/)]). 



Note that, as said in Appendix D, Eq.(4.15) holds only if the multipoles are evaluated 
with respect to world-lines w^{t) — z^{t, f]{T)) with f]{T) — ff — const., namely with respect 
to one of the integral hues of the vector field zif{T, a) dfj,. 

For an isolated system described by the multipoles appearing in Eq.(4.15) [this is not 
true for those in Eq.(4.4)] the equations df^T^^^O [see Eqs.(D4) and (D7)] imply no more 
than the following Papapetrou- Dixon- Souriau equations of motion [27, 68, 77, 78] for the 
total momentum P!^(Ts) — e'^{u{ps))qT^ (Tg) ~ and the spin tensor St^(Ts)[fj — 0] 



dT. 



0, 



dT, 



2P|^(T,)w^J(p,) = 2Klel^ {u{ps)X\ps) ^ 0, 



or ^-0, ^-0, ^^0. (4.16) 
dT, ' dT, ' dT, ^ ^ 
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V. DIXON'S MULTIPOLES AND RELEVANT CENTROIDS FOR CLOSED AND 
OPEN SYSTEMS OF INTERACTING RELATIVISTIC PARTICLES. 



In this Section we present new applications of the multipolar expansion to interacting 
systems of particles and fields. We first deal with the case of an isolated system of positive- 
energy relativistic particles with mutual action-at-a-distance interaction (see Section VIII 
of Ref. [2] and Section VI of Rcf.[5]); then wc deal with the case of an open particle sub- 
system of an isolated system consisting of charged positive-energy relativistic particles 
(with Grassmann-valued electric charges to regularize the Coulomb self-energies) plus the 
electro- magnetic field [5]. 



A. An isolated System of Positive-Energy Particles with Action-at-a-Distance In- 
teractions. 

As said in Section VIII of Rcf. [2], in the rest-frame instant form the most general expres- 
sion of the internal energy for an isolated system of N positive-energy particles with mutual 
action-at-a-distance interactions is 

M = Y^^/^^if^U~V(^~~^^ + V, (5.1) 

i 

— * 

where all the potentials Ui, Vi, V are functions of Ki- Kj, Ifji—ffjl, i<-k ■ {Vi — Vj) ■ the other 
hand, as shown at the end of Section II, in the free case we have 



•^ifree) = + N Ibi T^qa ' T^qb- (5.2) 



Since the 3-centers R+ and q+ became interaction dependent, in the interacting case we 
do not know the final canonical basis /?+, pqa-, Tiqa explicitly. For an isolated system, 
however, we have M = \J M? -|- ~ M. with M. independent of ({M, /?_!_} = in the 
internal Poincare' algebra). This suggests that also in the interacting case the same result 
should hold true. Indeed, by its definition, the Gartenhaus-Schwartz transformation gives 
Pqa ~ Pai T^qa ~ also in prcscncc of interactions, so that we get 



= MU+=o = V ^? + + (i^i - Vi)'' + V, (5.3) 

i 

where the potentials Ui, Vi, V are now functions of vf^a • 7?g6, T^qa ■ Pqb, Pqa ■ Pqb- 
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A relevant example of this type of isolated system has been studied in Ref.[5] starting 
from the isolated system of charged positive-energy particles (with Grassmann-valued 
electric charges Qi = 9* 9, = 0, QiQj = Qj Qi ^ for i ^ j) plus the electro- magnetic 
field. After a Shanmugadhasan canonical transformation, this system can be expressed 
only in terms of transverse Dirac obscrvables corresponding to a radiation gauge for the 
electro-magnetic field. The expression of the energy-momentum tensor in this gauge will be 
shown in the next Subsection. In the semi-classical approximation of Ref.[5], the electro- 
magnetic degrees of freedom are re-expressed in terms of the particle variables by means of 
the Lienard-Wiechert solution in the framework of the rest-frame instant form. In this way 
it has been possible to derive the exact semi-classical relativistic form of the action-at-a- 
distance Darwin potential in the reduced phase space of the particles. Note that this form is 
independent of the choice of the Green function in the Lienard-Wiechert solution. In Ref . [5] 
the associated energy-momentum tensor for the case N — 2 [Eqs.(6.48)] is also given. The 
internal energy is M = ^/M^T^ ^M^ E?=i V'^i +^^ + [1 + • |)] where 
V is given in Eqs.(6.34), (6.35) [in Eqs. (6.36), (6.37) for mi = m2]. The internal boost K, 
[Eq.(6.46)] allows the determination of the 3-center of energy i?+ — in the 

present interacting case. 

The knowledge of the energy- momentum tensor T"^^ (r, a) and of ~ g+ allows to 
apply our formalism to find the barycentric multipoles of this interacting case. It turns out 
that, in the gauge Rj^ f=:i g+ ~ y+ ~ 0, all the formal properties studied in the previous 
Section (like the coincidence of all the relevant centroids) are reproduced in presence of 
mutual action-at-a-distance interactions. 



B. Open Subsystem of the Isolated System of N Positive-Energy Particles with 
Grassmann-valued Electric Charge plus the Electromagnetic Field. 

Let us now consider an open sub-system of the isolated system of A^ charged positive- 
energy particles plus the electro-magnetic field in the radiation gauge. The energy- 
momentum tensor and the Hamilton equations on the Wigner hyper-plane of the isolated 
system are, respectively, = Ylii ^i') — means evaluated on the equations of motion; to 
avoid degenerations we assume that all the masses are different] 
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^ I 



i=l 



1 



N 



i=l 



i=l 

TV 



(9 1 

•i=i 



l..Ar 



N 



T^^{T,a) = J2S%a-ff,{T)M{T)-Q,Al{T,f},{T))] + 



i=l 



+ [(7fx + 5^g,^5^(a-7/,(r))) x^](T,a), 

i=l 



i=l 



1 ^ 1 2 

-<5-[(7f^ + ^g,-53(a-rf,(r))) + ^ V 

N ^ N ^ 

i=l 1=1 



i=l 



(5.4) 
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A±_r{T,a) = -7r_Lr(T, ct), 



.^(r) + j dM^..B]ir,B)^Hrest-framecondmon). (5.5) 



Let us note that in this reduced phase space there are only either particle-field interactions 
or action-at-a-distance 2-body interactions. 

The particle world-lines are (r) = x'^-\-u'^{ps) T-|-e(f(w(ps)) 77[(t), while their 4-momenta 
are (r) = ^mf + - ^(r, 77,)]^ ^^'^(p.) + e(^(^i(p.)) K - ^K^, ^«)]- 

The generators of the internal Poincare group are 



1=1 

2 ^47r I r/,(r)-r/,(r) I J 2^ ^ ' ^' 

P(^int) = fi+{T) + J d^a[n± X B]{T,a) ^ 0, 



i=l 

N 



j=i 

^ AT I..N „ 

+ 2 5Z S / ""^^ ~ "^'^^^^ ' '^^'^ ~ "^'^^^^ + 

i=i j^i 

+ Q, J d'a7rl{T,a)c{a- Vi{r))] - ^ j d' a {ttI + B'){t, a), (5.6) 
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with c{fii - fjj) = l/(47r|77,- - fji\) [Ac(a) = 5^{a), A = -9^, c{a) = dc{a) = a/(47r 

Note that 'P^^j) = Q'^^ and 'P('j„() = ^^^^ are the mass and momentum monopoles, respec- 
tively. 

For the sake of simphcity let us consider the sub-system formed by the two particles of 

mass nil and m2. Our considerations may be extended to any cluster of particles both in 
this case and in the case discussed in the previous Subsection. This sub-system is open: 
besides their mutual interaction the two particles have Coulomb interaction with the other 
N — 2 particles and feel the transverse electric and magnetic fields. 

By using the multipoles we will select a set of effective parameters (mass, 3-center of 
motion, 3-momentum, spin) describing the two-particle cluster as a global entity subject to 
external forces in the global rest-frame instant form. This was the original motivation of the 
multipolar expansion in general relativity: actually replacing an extended object (an open 
system due to the presence of the gravitational field) with a set of multipoles concentrated 
on a center of motion. In the rest-frame instant form it is possible to show that there is no 
preferred centroid for an open system, namely different centers of motion may be selected 
according to different conventions unlike the case of isolated systems where, in the rest frame 

~ 0, all these conventions identify the same centroid. We will see, however, that there is 
a choice which seems preferable due to its properties. 

Given the energy- momentum tensor T^^(r, a) (5.4) of the isolated system, it would seem 
natural to define the energy-momentum tensor (r, a) of an open sub-system composed by 

a cluster of n < N particles as the sum of all the terms in Eq.(5.4) containing a dependence 
on the variables rfj, k,, of the particles of the cluster. Besides kinetic terms, this tensor 
would contain internal mutual interactions as well as external interactions of the cluster 
particles with the environment composed by the other N — n particles and by the transverse 
electro- magnetic field. There is an ambiguity, however: why attributing just to the cluster 
all the external interactions with the other N — n particles (no such ambiguity exists for the 
interaction with the electro-magnetic field)? Since we have 2-body interactions, it seems 
more reasonable to attribute only half of these external interactions to the cluster and 
consider the other half as a property of the remaining N — n particles. Let us remark that 
according to the first choice, if we consider two clusters composed by two non-overlapping 
sets of Ui and n2 particles respectively, we would get ^cf?fi+n2) -^cfjfi) + -^cfn^)' since the 
mutual Coulomb interactions between the two clusters are present in both and T^^^y 

Instead according to the second choice we would get T4^ , „ n = T45n -|- T^\. Since this 
property is important for studying the mutual relative motion of two clusters in actual 
cases, we will adopt the convention that the energy-momentum tensor of an particle cluster 
contains only half of the external interaction with the other N — n particles. 

Let us remark that, in the case of fc-body forces, this convention should be replaced 
by the following rule: i) for each particle mj of the cluster and each /c-body term in the 
energy-momentum tensor involving this particle, we write k = hi -\- {k — hi), where hi is 
the number of particles of the cluster participating to this particular A;-body interaction; ii) 
then only the fraction hi/k of this particular k-hody interaction term containing rrii has to 
be attributed to the cluster. 

Let us consider the cluster composed by the two particles with mass mi and 777.2. 
The knowledge of T^^ = T^S on the Wigner hyper-plane of the global rest-frame in- 
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stant form allows to find the following 10 non conserved charges [due to = we have 



2 

Ql Q2 1 Qk _ 

+ 47r |r/i(r) - rf2(r)P 2 ^ 47r |77,(r) - 77,(r)P ' 

{y"dVTnT,a)} = ^i(T)+^2(T), 

= ffi{T) X Ki(r) +r;2(r) x K2(t), 

iC^^ - j d^aaT^^{T,a) = 
2 



^i(^) V"^i + ['^ilT-) - QiMr,f]i{r)W - 

i=l 

^Qi d^aTf±{T, a) c{a - f)i{T)) - 
Q1Q2 J d^aac{a-f]i{T))-c{a-f]2{T))- 
\Y Qi Y Qk J d^aac{a - f]i{T)) ■ c{a - ffk{r)) 



i=l k^l,2 



K,c{int) = - r)i{T) \l mj + kI - Qi Q2 / d^a ac{a - f]i{T)) ■ c{a - f]2{T))l5.7) 

i=i 

which do not satisfy the algebra of an internal Poincare' group due to the openess of the 
system. Since we work in an instant form of dynamics only the cluster internal energy and 
boosts depend on the (internal and external) interactions. Again Mc = ql'^ and = q^J 
are the mass and momentum monopoles of the cluster. 



Another needed quantity is the momentum dipole 
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2 2 „ 

i=i 1=1 

2 

-2Y,Q, I ciVc(a-7f,(r))[aMi + aMl](r,a), 



The time variation of the 10 charges (5.7) can be evaluated by using the equations of 
motion (5.5) 
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i=l 



+ 



0. 



— * 



2 



"if + «f 



47r|7;i(r) -7ffe(r)|3 



X >lx(T,?7i(T)) +?7i(T) X 



^±(r,77,(r))- 



- E^^ 



r/i(r) X rffc(r) 
47r |r;i(r) - ^(r)^ 



2 

- E 

i=l 
2 

2 



7r±(T, 77i(T)) + E <5fe c(77j(T) - 77fe(T))J + 



ciViir) - ffkir)) - / rfV7r_L(r,a) 



Ki(r) ■ c(a - rfi(r)) 



+ 



1=1 fc^i 
Q1Q2 / (iV<T 



Ki(r) 



+ cl 



c{a - f]i{T)) ■ c{a - f]2{T)) + 



- 9 E E / c^'cr<? 



i=l fc^l,2 



+ 



c(a-77.(r))- [(^ f^^^^ ■a)ci:a-77,(r))]). (5.9) 



Let us remark that, if we have two clusters of n\ and particles respectively, our 
definition of cluster energy-momentum tensor implies 



— * — * — * 

^c(ni+n2) ^c(ni) "I" ^0(112)7 
— * — * — * 

ni+n2) — c7c(ni) ~l~ c7c(n2)) 

^c(ni+n2) — ^c(ni) + ^c(n2)- (5.10) 
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The main problem is the determination of an effective center of motion Q{r) with world- 
hne w^{r) = + it^(Ps) r + e!^{u{ps)) Q{r) in the gauge Tg = r, q+ ^ R+ ^ y+ = of the 

/ ^ 

isolated system. The unit 4-velocity of this center of motion is u'^.{t) = w^{t)/\/ 1 — C^{t) 
with w^{t) = u^{ps) + e^(M(Ps)) Cc{t)- By using 5 z^{r, a) = e'^{u{ps)) (cr'' — C{t)) we can 
define the multipoles of the cluster with respect to the world-hne u'^(r) 

qr.-ruAB^^^ = J d'a [<7- - Q{T)]..[a^- - C"(r)] T^^^ir, a). (5.11) 

The mass and momentum monopoles and the mass, momentum and spin dipoles are 
respectively 



= Pc"(^)-C(r)P«, 

Sr = [e'^Hps)) W'ip,) - e';{u{ps)) u>^{ps)] q7^ + 6^(«(p.)) el{u{ps)) (g^ - C^) 
= [e^(«(p.)) «''(p.) - e;: («(!>,)) «np,)] M, (7?^ - C^) + 



m 



ciPs) = --^r MPs) = -e^rHPs)) ?r ■ (5.12) 



Let us now consider the following possible definitions of effective centers of motion (many 
other possibilities exist) 

— * — * — * 

1) Center of energy as center of motion, C,c{e){t) = Rc{t), where Rc{t) is a 3-center of 
energy for the cluster built by means of the standard definition 

^ = -|. (5.13) 

It is determined by the requirement that either the mass dipole vanishes, q^"^ = or the 
mass dipole moment with respect to u^{ps) vanishes, '^^i^s) ~ ^' 

The center of energy seems to be the only center of motion enjoying the simple compo- 
sition rule 

w Mc(„i) ^c(ni) + Mc(n2) Rc{n2) , - ^ 

Rc{ni+n2) = 77 ■ (5-14j 

c(ni+n2) 

— * — * 

The constitutive relation between Vc and Rc{t), see Eq.(4.8), is 
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= 



dq, 



TTT 



dr 



-ICl- McRl- M,Rl, 



Vc = M,R, + M,R,- 

2 



1=1 
2 



7r±(T, f]i(T)) + ^Qk c{fii(T) - ffkir)) 



+ 



'?fc(T) 



Ki{T) ■ c{a - f]i{T))- 



%=\ k^i ^y'^k^'^kV^) 



+ 



= <9i(92 / d^aa 



2 



d^c{B - t]x{t)) ■ c{a - if2{T)) + 



i=l k=/=l,2 



J^Q^ E Q'^ d'aa(\(^^^=-d)c{a-fi,{r))]-c{a-fik{^^^ 



+ cia - f)AT)) ■ \( , ""^^^^ ■ d] c(a - ^fr))!). 



(5.15) 



Prom Eq.(4.10) the associated cluster spin tensor is 



dc J 



= e(^(«(p,)) 6:(«(p,)) 6^™ \j: - {R, X P,)^ 



(5.16) 



2) Pirani centroid Cc{p){t) (is center of motion. It is determined by the requirement 
that the mass dipole moment with respect to 4- velocity w^(t) vanishes (it involves the 
anti-symmetric part of p"*") 



c{wo) = --ST ^ci. = 0, =^ Cc(p) • Cc(p) = Cc(p) • Rc, 



Cc(P)(t) 



M,-P,-Ce(p)(T) 



M,i?,-it:,-C(p)(T)P,-Cc(P)(T) X J-J. (5.17) 
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Therefore this centroid is imphcitly defined as the solution of these three coupled first 
order ordinary differential equations. 



3) Tulczyjew centroid Cj{t) (t) as center of motion. If we define the cluster 4- momentum 

= Mcui'{ps)+Vl ei^{u{ps)) [P^ = Ml-Vl M% its definition is the requirement that 
the mass dipole moment with respect to P^ vanishes (it involves the anti-symmetric part of 

p7) 



m 



c{Pc) 



Pcu = 0, ^Vc- Cc(T) = Vc ■ Rc, 



Cc{T){t) 



Ml - VI 



Ml Rc-Vc- RcVc -VcXJc 



(5.18) 



Let us show that this centroid satisfies the free particle relation as constitutive relation 

Vc = McCciT), 



P!f = Mr 



ri 1 

1c(T) 



M, 



M2 - P2 



V' Rc + Vc- RcVc + Vc X Jc 



+e(f (i.(p,)) el{u{ps)) e™'' \j: - (Cc(t) x Vcf 



(5.19) 



— # — * 

If we use Eq.(5.17) to find a Pirani centroid such that C^. = Vc/Mc, it turns out that the 
condition (5.17) becomes Eq.(5.18) and this implies Eq.(5.19). 

The equations of motion 



Mc{r) Cc(T)(r) = Vc{t) - Mc{t) CciT)ir), 



(5.20) 



contain both internal and external forces. Notwithstanding the nice properties (5.19) and 
(5.20) of the Tulczyjew centroid, this effective center of motion suffers the drawback of not 
satisfying a simple composition property. The relation among the Tulczyjew centroids of 
clusters with ni, n2 and ni + 77,2 particles respectively is much more complicated of the 
composition (5.14) of the centers of energy. 

All the previous centroids coincide for an isolated system in the rest-frame instant form 
with Vc = K'+ ^ in the gauge g+ ~ R+ a; y+ a; 0. 
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4) The Corinaldesi-Papapetrou centroid with respect to a time-like observer with 4-velocity 

Cilcp)i'^) center of motion. 

m^) = -Srvu = 0. (5.21) 

Clearly these centroids are unrelated to the previous ones being dependent on the choice 
of an arbitrary observer. 

5) The Pryce center of spin or classical canonical Newton- Wigner centroid Cc{nw)- 
It defined as the solution of the differential equations implied by the requirement 
{Q{NwyCc{NW)} = 0, {C(w)'^c} = Let us remark that, being in an instant form 
of dynamics, we have {VI, V^} — also for an open system. 

The two effective centers of motion which look more useful for applications seems to be 

— * — * — * 

the center of energy Cc{e){t) and Tulczyjew's centroid Cc{t){t), with Cc{e){t) preferred for 
the study of the mutual motion of clusters due to Eq.(5.14). 

Therefore, in the spirit of the multipolar expansion, our two-body cluster may be de- 
scribed by an effective non-conserved internal energy (or mass) Mc{t), by the world-line 
= -I- u^{ps) r -\- €lf{u{ps)) Q(EorT)(''') ^ssociated with the effective center of motion 

Cc(EorT){T) and by the effective 3- momentum Vc{t), with Cc{EorT){T) and Pc(t) forming a 
non-canonical basis for the collective variables of the cluster. A non-canonical effective spin 
for the cluster in the 1) and 3) cases is defined by 



a) case of the center of energy, 



Sc{E){r) = Jc(t) - Rc{t) X Vc{t), 



dCc[E){T) _ dXjr) dR^jr) ^ ^ dV^jr) 



b) case of the Tulczyjew centroid, 

5c(t)(t) = Jc{t) - Cc(r)(T) X V,{t) = 

^ M'At) S,^e){t) - Voir) ■ Xjr) Vc{t) 

Since our cluster contains only two particles, this pole-dipole description concentrated 
on the world- line w^{t) is equivalent to the original description in terms of the canonical 
variables rfiir), ^^(t) (all the higher multipoles are not independent quantities in this case). 
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Let us see whether it is possible to replace the description of the two body system as an 
effective pole-dipole system with a description as an effective extended two-body system by 
introducing two non-canonical relative variables Pc{EorT){T), 7^c{EorT){T) with the following 
definitions 



?7l = U{Em-T) + 2 PciEorT), U{EorT) ^ [Vl + V2), 

^ def -> 1^ ^ _ ^ ^ 

V2 — U{EorT) — -^PciEorT), Pc(EorT) — Vl ~ V2: 

-> def 1 ^ ^ ^ ^ ^ 

«1 = -l-^c + T^ciEorT), Pc = K.1 + 1^2, 

K'2 = TTc(EorT), T^c{E orT) ^ - 1^2) , 

— * — * — * 

Jc = 77*1 X 77*2 X = Cc{EorT) X Vc + Pc{EorT) X T^c{EorT), 

=^ Sc{EorT) — Pc{EorT) X T^c{EorT)- (5.23) 



Even if suggested by a canonical transformation, it is not a canonical transformation and it 
only exists because we are working in an instant form of dynamics in which both Vc and Jc 
do not depend on the interactions. 

Note that we know everything about this new basis except for the unit vector 

Pc{EorT)/\Pc(EorT)\ and the momentum T(c{EorT)- The relevant lacking information can be 
extracted from the symmetrized momentum dipole + p^^ . which is a known effective 
quantity due to Eq.(5.9) having the following expression in terms of the variables (5.23) 



p™+pr + 25]gi / dVc(a-77,(r))rAl + 9^Al](T,a) = 

i=l 

2 

^ ^ ^^"i '^i + '^i ) ^ Q(EorT) '^c + Q{EorT) + 
1=1 

"T Pc{E or T) i^c^E or T) "r Pc{E or T) c(E or T) ■ 

(5.24) 

A strategy for getting this information is to construct a spin frame, which, following 
Ref.[2] for the N = 2 case, is defined by Sc{EorT) = Sc(EorT)/\Sc{EorT)\, T^dEorT), Vc(EorT) = 

T^c{EorT) xSc{EorT)-, '^^i^Sc{EorT)-T^c{EorT) = ^c{EorT) ^ ^^EorT) ^ ~^c(EorT) ^ Then 

we get the following decomposition {{Sc(EorT) = \Sc(EorT)\) 
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Pc{EorT) — Pc(EorT)'^c(EorT), Pc(£ or T) — l^fl ~ ''?2 1 , 

~ ^ Sc{EorT) . ^ Pc{EorT) 

T^c{EorT) — T^c{EorT) l<-c{EorT) Vci^EorT), T^c{EorT) — T^c{EorT) ' , 

Pc{EorT) Pc{EorT) 

(5.25) 

where Pc(EorT) is just the relative variable appearing in the Coulomb potential. Eqs.(5.25) 

show that only the three variables Tic{EorT) and HciEorT) — Pc{EorT)/ Pc{EorT) are still un- 
known. Then from Eqs.(5.23) and (5.24) we get 



PI{E or T) ^c(E or T) + Pc(E or T) ^1{E orT) — 2 Pc{E or T) TTc(E or T) T^l(E or T) "^"(B or T) ~ 



— Sc{E or T) ( T^'c^E or T) '^c(E or T) 



'^c{EorT) '^c{EorT) 



— Pc" +PT ~ {Cc{EorT) orT) 

v:) + 

iQi I dVc(a-77;(T))[9'-Al + 9"A;](T,a) 
1 



2 



c{E or T) ) 



K{EorT)^c{EorT) = 0- (5.26) 

But these are three independent equations for fciEorT) and for the two degrees of freedom 
in the unit vector H-ciEorT) in terms of the known quantities F'^{EorT)i ^c(EorT), Pc[EorT) — 
\fji — fj2\- For instance we get iTc{EorT) = (^^ ^c(EorT))/'^ Pc{EorT)'- due to the transversality 
of the vector potential TTc(EorT) does not depend on it. In conclusion the external electro- 
megnetic potential A± enters only in the determination of the axis 'R'c(EorT) of the spin 
frame. 

This completes the construction of the effective relative variables and of the effective spin 
frame using the extra input of the 3-momentum dipole. In this way we get a description 
of the two-body cluster as an effective two-body system instead of a pole-dipole system. 
However, the weak point of this description of the open system as an extended object is 
that, whatever definition of effective center of motion one uses, the symmetrized momentum 
dipole + p^^ does not depend only on the cluster properties but also on the external 
electro- magnetic transverse vector potential at the particle positions, as shown by Eq.(5.8). 
As a consequence the spin frame, or equivalently the 3 Euler angles associated with the 
internal spin, depends upon the external fields. 

If we accept this drawback, it is reasonable that, by taking into account higher multipoles, 
it is possible to give a description of a cluster of n > 3 particles in terms of as many effective 
n-body systems as effective dynamical body frames following the scheme of Ref.[2]. 

This would open the possibility to have effective descriptions of two clusters of Ui and 
n2 particles, respectively, and to compare it with the effective description of the cluster 
composed by the same rii + n2 particles to find the relation between the three centers of 
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motion of the (ni +^2)-, ni- and 77,2- clusters and the relative motion of the two rii- and ^2- 
clusters. To this end the use of the center of energy as center of motion seems unavoidable 
due to the simple composition law (5.14). Whatever choice we adopt, however, it turns 
out that the relative motion of the two clusters depends on the external fields besides the 
effective parameters of the clusters. 

These techniques can be extended to relativistic perfect fluids, if described in the rest- 
frame instant form as done in Refs. [1]. Moreover, they are needed for the determination 
of the post-Minkowskian approximation to the quadrupole formula for the emission of grav- 
itational waves (re-summation of the post-Newtonian approximations) in the background- 
independent Hamiltonian linearization of tetrad gravity [79] plus a perfect fluid [80] . 
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VI. CONCLUSIONS. 



A relativistic description of open systems like binary stars embedded in the gravitational 
field would be an important achievement nowadays in view of the construction of templates 
for the gravitational radiation. Even by approximating such description by means of a 
multipolar expansion in a way suitable for doing actual calculations, either analytical or 
numerical, a big amount of kinematical technical preliminaries is needed anyway. With this 
in view, we had in mind to develop methods which could be useful in general relativity with 
relativistic perfect fluids as matter, where single or binary stars could be described by open 
fluid subsystems of the isolated system formed by the gravitational fleld plus the fluid in the 
rest-frame instant form of either metric or tetrad gravity [79, 80]. 

To pursue our program, in the present paper we have first of all completed the study 
of the relativistic kinematics of the system of N free scalar positive-energy particles in the 
rest-frame instant form of dynamics on Wigner hyper- planes, initiated in Ref.[2]. 

Then, we have evaluated the energy momentum tensor of the system on the Wigner 
hyperplane and then determined Dixon's multipoles for the N-body problem with respect 
to the internal 3-center of mass located at the origin of the Wigner hyperplane [81]. For 
an isolated system most of the existing definitions of a collective centroid identify a unique 
world-line, associated with the internal canonical 3-ccntcr of mass. In the rest-frame instant 
form these multipoles are Cartesian (Wigner- covariant) Euclidean tensors. While the study 
of the monopole and dipole moments in the rest frame gives information on the mass, the 
spin and the internal center of mass, the quadrupole moment provides the only (though not 
unique) way of introducing the concept of barycentric tensor of inertia for extended systems 
in special relativity. 

By exploiting the canonical spin bases of Refs.[2, 21], after the elimination of the internal 
3-center of mass {q+ = k+ = 0), the Cartesian multipoles q!^ ' '^"^^ can be expressed in 
terms of 6 orientational variables (the spin vector and the three Euler angles identifying the 
dynamical body frame) and of 6A^ — 6 {rotational scalar) shape variables, i.e. in terms of the 
canonical pairs of a canonical spin basis. 

Having completed the discussion of the isolated system of positive energy free scalar 

particles the previous formalism has been applied to an isolated system of positive-energy 
particles with mutual action- at- a- distance interactions. Here again we find a unique world- 
line describing the collective motion of the system. 

On the other hand, in the case of an open n < N particle subsystem of an isolated 
system consisting of A^ charged positive-energy particles plus the electro-magnetic field a 
more complex description appears. In the rest frame of the isolated system a suitable 
definition of the energy-momentum tensor of the open subsystem allows to define its effective 
mass, 3-momentum and angular momentum. However, unlike the case of isolated systems, 
each centroid putatively describing the collective centers of motion, gives rise to a different 
world-line. Starting from the evaluation of the rest-frame Dixon multipoles of the energy- 
momentum tensor of the open subsystem with respect to various centroids we are given 
therefore many candidates for an effective center of motion and for an effective intrinsic 
spin. Two centroids (the center of energy and Tulczyjew centroid) seems to be preferred 
because of their specific properties. In the case n = 2 it is possible to replace the pole- 
dipole description of the 2-particle cluster with a description of the cluster as an extended 
system (whose effective spin frame can be evaluated) at the price of introducing an explicit 
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dependence on the action of the external electro-magnetic field upon the cluster. 

Finally, by comparing the effective parameters of an open cluster of rii +n2 particles with 
the effective parameters of the two clusters with ni and n2 particles, it is shown that only 
the effective center of energy can in fact play the role of a useful center of motion. 

The kinematical concepts we have defined for closed and open N-body systems are enough 
for the treatment of relativistic continua like relativistic fiuids. In [1] a preliminary extension 
to closed relativistic fiuids is given. 
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APPENDIX A: NON-RELATIVISTIC MULTIPOLAR EXPANSIONS FOR N 
FREE PARTICLES. 



In the review paper of Ref.[27] it can be found a study of the Newtonian multipolar ex- 
pansions for a continuum iscntropic distribution of matter characterized by a mass density 
p{t,a), a velocity field If^itja), and a stress tensor (T^^{t,a), with p{t,(x)U{t,(T) the mo- 
mentum density. In case the system is isolated, the only dynamical equations are the mass 
conservation and the continuum equations of motion, respectively 



dp{t,a) dp{t,a)U'{t,a) _ 

dt da^ 
dp{t, a)U''{t, a) _ d[pU^U' - a^']{t, a) ^ 

dt da' ~ 



0. (Al) 



We can adapt this description to an isolated system of N particles in the following way. 
The mass density 



iV 



p{t,a) = Y,^^^\^-V^{i)). (A2) 



1=1 



satisfies 



while the momentum density [82] is 



N 



p{t,a)U^{t,a) = 5^m,7/,(i)(53(a-7f,(t)), (A4) 



i=l 



The associated constant of motion is the total mass m = 

If we define a function ({(T,ffi) concentrated in the N points ffi, i=l,..,N, such that 
rfj) = for (T 7^ f/j and CiVi^Vj) = ^ij [83], the velocity field associated to N parti- 



cles becomes 

N 



i=l 

The continuum equations of motion are replaced by 



8 8 ^ ^ 



dt 



i=l 



i=l 



dg djpU'-W -a'']{t,a) 
da' 



(A6) 
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For a system of free particles we have f]^{t) = so that a'^^{t, a) — 0. If there are inter- 
particle interactions, they will determine the effective stress tensor. 

Let us consider an arbitrary point rj{t). The multipole moments of the mass density p 

— * 

and momentum density pU and of the stress-like density pU^U^ with respect to the point 
ff{t) are defined by setting {N > 0) 



m''-'"[f]{t)] = J ciV[cr"i - r?''^(t)]...[cr"" - r/''"(t)]p(r, a) = 

N 
i=l 

N 

n = m[f]{t)] = m = ^^mj, 

1=1 

N 

= J2^^r)Kt)[r|l^it) - ri^^{t)]...[rjl''{t) - v'^it)], 

i=l 

N N 

n = p'^ivit)] = ^Mt) = XI 4 = ~ 0, 

i=l i=l 

p^^--r.rs^fj{t)] = j - r]'\t)]...[a^- - 77^"(i)]p(i, BW{t, a)U%t, a) = 

N 

= Y.''^^mm[v?{t) - r,^^{t)]..M-{t) - (A7) 



i=l 



The mass monopole is the conserved mass, while the momentum monopole is the total 
3-momentum, vanishing in the rest frame. 

If the mass dipole vanishes, the point f]{t) is the center of mass: 



N 



i=l 



The time derivative of the mass dipole is 

'^"'^]'^^^^^ = p''[fj{t)] - mf}\t) = k\- mff{t). (A9) 

When f]{t) = (jnr, from the vanishing of this time derivative we get the momentum-velocity 
relation for the center of mass 

P^[Qnr] = = mg!j_ [fti in the rest frame]. (AlO) 
The mass quadrupole is 
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N 
i=l 

SO that the barycentric mass quadrupole and tensor of inertia are, respectively 



N 



i=l 

u 

i=l 

1...N-1 



a,b 

N-1 

m''[qnr] = S^' Yl ^-bPa ' Pb - 1''%^. (Al2) 
a,6=l 

The antisymmetric part of the barycentric momentum dipole gives rise to the spin vector 
in the following way 

N N 



i=l i=l 
N-1 



a=l 



The multipolar expansions of the mass and momentum densities around the point fflt) 
are 



p(t,a) = > S\a-f^{t)), 



n=0 



°° ^'•i....rnrr^ Qn 



71=0 

Finally, for the barycentric multipolar expansions we have 
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pit, a) = m6%a - q^r) - l{r'[Qnr] - Yl ' 5-) + 



+ E 



^r-i....r„J-^ 



n=3 



d_ 

J 



pri....r„rj- 



n=2 



(A15) 
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APPENDIX B: SYMMETRIC TRACE-FREE TENSORS. 



In the applications to gravitational radiation, irreducible symmetric trace-free Cartesian 
tensors (STF tensors) [76, 84, 85, 86] are needed instead of Cartesian tensors. While 
a Cartesian multipole tensor of rank I (like the rest-frame Dixon multipoles) on E? has 3' 
components, |(/ + l)(/ + 2) of which are in general independent, a spherical multipole moment 
of order I has only 21 + 1 independent components. Even if spherical multipole moments 
are preferred in calculations of molecular interactions, spherical harmonics have various 
disadvantages in numerical calculations: for analytical and numerical calculations Cartesian 
moments are often more convenient (see for instance Rcf. [87] for the case of the electrostatic 
potential). It is therefore preferable using the irreducible Cartesian STF tensors [88, 89] 
(having 2/ + 1 independent components if of rank /), which are obtained by using Cartesian 
spherical (or solid) harmonic tensors in place of spherical harmonics. 

Given an Euclidean tensor A].^ }.^ on i?^, one defines the completely symmetrized tensor 
Sk^,,ki = = ji XItt ^fc,r(i) - '^V(i)- Then, the associated STF tensor is obtained by 

removing all traces {[1/2] = largest integer < 1/2) 



[1/2] 

ASTF) _ A A Q 

^fei...fe/ ~ / ^n0{kik2---^k2n-ik2n'^k2„+i--ki)iiii...j„jrn 

n=0 



n = ( ir ^'(2/-2n-l)!! 
" - ^ ' (/-2n)!(2/-l)!!(2n)!!- ^ ^ 



For instance {TabcY'^^ = T^abc) - | SabT(^nc) + SacT^m) + SbcT^an)^ ■ 
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APPENDIX C: THE GARTEN AHUS-SCHWARTZ TRANSFORMATION. 

In Ref.[2] we defined canonical internal relative variables with respect to the internal 
3-center of mass by exploiting a Gartenhaus-Schwartz canonical transformation. The 
canonical generator of this transformation is 

G = q+- K+, (CI) 

so that the finite transformation, depending on a parameter a, on a generic function F on 
the phase space is 

F{a) = F+ da {F{a),G{a)}. (C2) 
Jo 

In particular we have 

lim K+{a) = 0, lim q+{a) — oo. (C3) 

a— »oo a— >oo 

As said in Section II, if we define the canonical transformation (2.19), then the quantities 

TTqa = lim 7f„(a), = lim Pa{a), (C4) 

a-^oo a-^oo 

are well defined and the transformation 

ffi, Ki K+, q+, Pqa, Tlqa, (C5) 

is a canonical transformation 

{g;, <} = <J = (C6) 

as said in Eq.(2.18). 

The quantities Pqa,7^qa are the searched internal relative variables: they describe the 
system after the gauge fixing q+ 0, k,+ ^ 0. We have also 

pqa^ pa, T^qa ~ TTa- (C7) 

Thanks to these results, we can calculate a function F independent of g+ on the phase 
space, under the constraint /«+ ~ 0, by simply performing the limit 

F U+«iO {pqa^T^qa) = 1™ F{a). (C8) 

a— >oo 

This method is applied in Section IV for calculating the multipoles after the gauge fixing 
g+ pa 0, 0. These multipoles depend on so that (see Ref.[2]) 

Af-l 



and on {fji — R+), so that 



lim Ki{a) = VNy^^-faiTfqa, (C9) 



a=l 
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m2 + ^2 



-EE - 7c.) A. . (CIO) 



Then using the (C9) and (C4) we have 



hm {fi,{a) - R^{a)) = E E ^(^^ - laj)Pqa J^ ^ =■ (Cll) 

The following notation is used to denote the hmits (C8) 

F —>-a^oo F |k_,_«o {PqaiT^qa)- (C12) 

For example Eqs. (C9) and (Cll) become 



7V-1 



a=l 



{fli-R+) ^a-.oo 2^2^VA^(7a»-7aj)Pqa „ / o , ^ ^ =" ^^'^^^ 

Efc V "^fe + N }2ab ^qa ' T^qblaklbk 



The closed form (2.20) - (2.22) of the canonical transformation (C5) was not given in 
Ref.[2], but it can derived from the following two equations of that paper [its Eq.(5.13) and 
(5.24)] 
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= na + —[{Msys - M)n+ ■ TTa - \K+\Ha] 



sys 



M2_ - 



TV 



AT 



N-l 



{rel)i 



a=l 



qa-i 



\ 



1..N-1 



ab 



i=l 
N 



N 



i=l 



L.W-l 



latybi^qa ' '^qb- 



ab 



def ^ / \ -> 
Pqa = Pa(oo) = p„ 



- ^^laj{lbi-lbj) 



i,j=l 6=1 

W N-l 



R+\'ilj{oo) ,M. 



+ 



sys 



-l)n+ 



• P6 ~ Pa- 



n+ ■ pb 



i,j=l b=l 



Msys Hj{oo)VTl 
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APPENDIX D: MORE ON DIXON'S MULTIPOLES. 



In this Appendix we shall consider multipoles with respect to the origin, i.e. with rf = 

[we use the notation ti^^'"^"^'^ (Tg) ttp^"''^"'^^(Ts, 0)] and we shall give some of their properties 
following Ref.[25]. The proofs of these results are identical to those given in Ref.[25] after 
the following kinematical modifications of Eqs.(3.6) and (3.7) of that paper: 

i) let W be the world-tube containing the compact support of the isolated system and 
w^{t) = z^^ij, fjij)) a time- like world- line inside it, with tangent vector w^{t) = z^{t, + 
z^{t, f]{T)) rf{T), used to evaluate the multipoles (see Section IV); 

ii) the Wigner hyper-planes T^wt of the rest-frame instant form are in general not orthog- 
onal to the world-line (differently from the hyper- surfaces S(s) of Rcf. [25]); 

iii) Eq.(3.6) is replaced with J d^z f\z) = J dr d^a a/ g{T, a) f{z{r, a)) — 
J dr d^'EnZl^{T,a) f{z{T,a)) (see before Eq.(2.1) for the notations); 



iv) since /^(r,a) = [Ip^A^tK^,^) = [J7,zJ,]{r,a 



we 



get J^d^f{z{T,a)) = J^^ d^a ^/-f{T, a) 1^{t, a) f{z{T, a)) 

j^^ d^a ^/'g{T\¥jzl{T,a) f{z{T,a)), so that we ha ve £ J^^ d^E,, f{z{T,a)) 
kr d"^^ Wair, a) zI{t, a ) f{z{r , a))] = f^^ d'^a Oa Wg{T, a) z^{t, a) f{z{T, a))] after 
having added J^^ d^a dr [\/g{r, a) z'^{t, a) f{z{T, a))] = 0; 

v) since ^\^zi]{T,a) = £^ [1, e^^^^ e,p,s §^ §^ ^]{r,a) = 0, we get that 
Eq.(3.7) is repla ced by f ci% /(z(r, a)) = ^ d^a z^{t, a) ^^^^^ = 

/^^ V7(^/.(r, a) <(r, a) = J^^ d^K <(r, a) '-I^; 

vi) as a consequence the results quoted in this Appendix hold if the multipoles are taken 
with respect to world-lines w^{t) — z^{T,f]) such that ?7(r) — if — const, (they are the 
integral lines of the vector field z!^{t, a) d^), so that w'^{t) — z!^{t, ff). 

As shown in Ref . [25] , if a field has a compact support W on the Wigner hyper-planes T,wt 
and if f{x) is a C°° complex- valued scalar function on Minkowski space-time with compact 
support [90], we have 



< T'^^ / > = j d^xT^''{x)f{x) = 

= j dT, j d'af{xs + 5xs)T^^^[xs{T,) + 5a;,(a)][0] = 

dT^j d'a I ^/(fc)e-'=-[-^(^^)+^-^(^)lr''^[a;,(r,) + 5a;,(a)][0] = 

n=0 

/ / 7C^/(^)^"*"^'^^^ E ^^K,...k,X--^-^''{Tsl (Dl) 
J J (27r)4 ^ n\ 



= / dT. 
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If, in particular, f{x) is analytic on W [25] [91], we get 



n=0 



T'^i^) = E ^Tdxfdx^r. f dTs5\x - x.(T,))t--— (T,). (D2) 

n=0 ■ 



For a N particle system this equation may be rewritten as Eq.(4.4). 
On the other hand, for non-analytic functions f{x) we have 



?i=0 

^ I J n=Ar+l 

and, as shown in Rcf. [25], from the knowledge of the moments t'^"'^'^^(Ts) for all n > N, we 
can get T'^'^{x) and, therefore, all the moments with n < N. 

In the case of N free particles and more in general for isolated systems, the Hamilton 
equations [92] for the multipoles (4.3), with p^^ •''"''(7;) ^^p^^ •''"^(7;, 0), imply 



,^ = 0, forn^O, 
dls 



If we define for n > 1 
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= e: 



Ml 



ri..r„Ar / 



1 



n+1 
n 



-w 



\Ps)q?-'^-^^{Ts) + e^MPsmp-'-'^^^iTs), 



and for n > 2 



(D5) 



n 



n + 2 



S^'-^'-''^'' {Ts) = 0, 



n — 1 



1 



+ 



+ -[u^{ps)e:{u{ps))+u''{ps)e^MPs))] 



n 



+ ^'sMPsMMPsWr-'-^'^'^T,) - 



n + 1 



+ 



(D6) 
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then Eqs.(D4) may be rewritten in the form 



1) n= 1 

N 

1^1 V + IS 
^^(T,) ^ 2p^{T,)vr^{p,)^2K%eY{u{ps))u''\ps)^^, 



2) n = 2 [identity f/-" ^4^^" + t^rf''^^ + t^rl^''^'] 



3) n > 3 

This allows to rewrite < T^*", / > in the form [2 5] 



< 



„ ft' 

n=2 



with 
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Irp 



iTs) 



r(Ml---Mn)(/i'^) 



i/il.../i„(UI' 



(Ts) 



n 
n 



n+1 



+ 



n 
n + 2 



n 



^Mn,/i ,1') 



Hps)) 



qr^l-rnAB^j.^^ _ 



n 



+ gj:^-'^"'^'^^(r,)]. (D9) 
For a N particle system, Eq.(D8)implies Eq.(4.15). 
Finally, a set of multipoles equivalent to the /^^ • is [93]: 



69 



for n >0 



^^ll■■■^lnl^l[pv]a]^rp ^ _ 



n+ 1 



!^]/Ul.../i„[p(T] 



{Ts) + 



(T]/il.../iri 
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(DIO) 



The Jl^^■■^^"^^'^P'^ are the Dixon 2"+^-po/e inertial moment tensors of the extended system: 
they (or equivalently the /^i- ''"^'^'g) determine the energy-momentum tensor together with 
the monopole and the spin dipole Sl^ . 

As shown in Section 5 of Ref . [25] , the equations of motion d^T'^'^ = do not imply equa- 
tions of motion for the multipoles j^i---^"'^^ (7^ > 2) or Jj^^■■■^^^^^^P'^ {n>Q), but only Eqs.(4.16) 
for the multipoles and Si^ [94]. Instead the multipoles f^^^■■■^^"^^'^ and p^i "^"^ have non 
trivial equations of motion. 

When all the multipoles jA^i- ^n/ii^po- ^^^^ negligible) one speaks of a pole-dipole 
system. 

On the Wigner hyperplane, the content of these 2'^^'^-pole inertial moment tensors is re- 
placed by the Euclidean Cartesian tensors q'^'"'^"'^'^ , q^'"'^"'^'^, q'^'"'^"^^. As shown in Appendix 
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B, we can decompose these Cartesian tensors in their irreducible STF (symmetric trace-free) 
parts (the STF tensors). 

Thus the multipolar expansion (4.4) may be rewritten as 



T^^[xf^^^{Ts) + e^Au{Ps)V] = T^''[wP{Ts) + 6f(«(p,))(a^ - r;^(T,))] 
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n 



+ e^^(«(p.))er,(«(p.))[g^,--^^^^(r„r/) - 



n 



leading to Eq.(4.15). 



(Dll) 



For open systems, subsystems of an isolated system like in Section V, we have 
d^T^^ = F^ _^ with an external force. As shown in Ref.[25] for the case in which 

^ i^Q^ jt^ =0) is the Lorentz force, in this case the multipolar expansion (4.15) 

is still valid, while the equations of motion (4.16) become {P^ and Tg are the conserved 
4-momentum and the rest-frame time of the global isolated system) 



2p^{Ts) u^\P) - j d'aa^ [e'^HP)) F^T^, a) - e^MP)) F^iUB^ 



dTs 

dsr{z)[ff=o] o 
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